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SECTION  I 


INTRODUCTION 


The  purpose  of  the  work  documented  in  this  report  is  to  improve  the 
accuracy  of  aerodynamic  hinge-moment  analysis.  When  addressing  this  subject., 
one  is  primarily  concerned  with  the  flow  over — and  the  pressure,  force, 
and  moment  experienced  by — the  nonplanar,  multiple  lifting  surfaces.  When 
the  control  surface  is  deflected,  the  pressure  distribution  on  the  lifting 
surface  is  significantly  altered  by  the  presence  of  the  discontinuous  kine- 
matic downwash  field.  Associated  with  the  change  of  geometry  are  such 
phenomena  as  boundary-layer  separation,  free  vortex  flow,  and  thick  wake. 

In  the  transonic  flow  regime,  the  shocks  lie  partway  back  on  the  lifting 
surfaces  and  interact  with  the  boundary  layer,  which  makes  the  problem  highly 
nonlinear.  Attention  to  these  phenomena  would  necessitate  a long-term  study 
and  a:  evaluation  of  each  individual  effect.  Indeed,  the  prediction  of  the 
hinge  moment  is  a complicated  subject.  In  view  of  the  current  status  of  the 
theoretical  development  and  the  numerical  schemes  in  various  aspects  of 
computational  aerodynamics,  one  might  want  to  place  reliance  for  such  informa- 
tion on  measurements  by  wind  tunnel  or  other  experimental  means  when  available. 

In  the  development  of  highly  maneuverable  aircraft,  the  accuracy  of  hinge- 
moment  prediction  is  important  with  respect  to  the  vehicle  stability  and 
controllability.  To  design  such  aircraft  certainly  calls  for  a more  efficient 
tool  using  the  product  of  theoretical  aerodynamics  instead  of  relying  solely 
on  measurements.  The  bulk  of  this  study  is  to  review  and  exemplify  several 
of  the  existing  methods  for  analyzing  the  aerodynamic  loadings  on  the  co-  ,.rol 
surfaces,  among  which  one  is  of  semiempirical  nature. 

The  first  part  of  this  work  (sec.  II)  is  concerned  with  reviewing  the 
formulation  and  coding  of  the  hinge-moment  subprogram  in  FLEXSTAB.  Aero- 
dynamic hinge-moment  derivatives  have  previously  been  computed  using  the 
FLEXSTAB  programs.  References  1 and  2 present  results  for  the  YF-16  and  B-52E 
aircraft,  respectively.  The  YF-16  analysis  was  conducted  by  the  AFFDL  and 
the  B-52E  analysis  was  completed  by  Boeing.  An  objective  of  che  present  study 
is  to  review  the  work  of  reference  1 and  to  improve  the  accuracy  of  hinge- 
moment  prediction. 


1 


Hinge  moment's  we^e  obtained  by  Boeing  for  the  YF-16  Government-furnished 
geometry  model.  The  effect  of  airplane  flexibility  on  hinge  moments  was 
investigated  using  AFFDL-supplied  NASTAP  data.  The  near-f ield/far-f ield 
option  was  applied  to  the  control  surfaces.  The  results  of  these  analyses 
were  compared  with  those  reported  in  reference  1.  To  investigate  the  effect 
of  paneling  on  hinge  moments,  the  YF-16  was  analyzed  with  various  paneling 
configurations  and  the  results  correlated  with  the  Government-furnished 
hinge-moment  data. 

The  second  part  of  the  work  (sec.  Ill)  is  to  review  and  discuss  the 
theoretical  and  numerical  characteristics  of  several  aerodvnamic  methods 
relating  to  the  load  prediction  on  the  lifting  surfaces.  These  methods 
include  FLEXSTAB,  Datcom,  RHO  4,  TEA  230,  and  che  PAN  AIR  pilot  code.  The 
discussions  are  focused  on  the  behavior  of  their  steady,  subsonic  solutions 
in  regions  of  singularity,  namely  at  the  control-surface  hinge  line  and  side 
edges,  and  at  planform  breaks  that  exhibit  discontinuities  in  surface  slope, 
planform  shape,  dihedral,  etc. 

Recommendations  are  made  in  sections  IV  and  V to  improve  the  aerodynamic 
methods  for  better  prediction  of  the  load  distributions  on  the  control  surfaces. 
A solution  procedure  that  combines  the  merits  of  panel  and  kernel  function 
methods  is  proposed  as  an  effective  scheme  to  predict  the  loading  on  a control 
surface  with  nonseparated  flow. 

Conclusions  are  presented  in  section  VI. 

SECTION  II 

FLEXSTAB  PROGRAM  REVIEW  AND  DEMONSTRATION 


2.1  HINGE-MOMENT  FORMULATION  AND  CODE 


The  hinge-moment  theory,  as  derived  by  Boeing,  was  reviewed  and  judged 
essentially  correct  in  the  context  of  the  aerodynamic  theory  as  used  in  FLEXSTAB 
(ref.  3,  vol.  I).  The  code  and  flow  charts  were  also  checked.  The  code  was 
found  to  describe  the  theory  properly;  however,  a few  minor  improvements  and 
corrections,  as  reported  in  monthly  progress  report  No.  2 (31  May  1977) 
submitted  to  AFFDL,  are  required  at  the  printout  level.  A large  number  of 
errors  were  found  in  the  hinge-moment  program  flow  charts  (ref.  3,  vol.  Ill) 
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2.2  HINGE-MOMENT  ANALYSIS  FOR  YF-16 


The  initial  aim  of  the  analysis  was  to  rerun  the  AFFDL  data  at  Boeing  to 
establish  a baseline  to  compare  the  AFFDL  results  with  the  Boeing  results. 
This  comparison  consisted  of  analyzing  the  basic  paneled  model  with  the 
residual  elastic  option  at  sea  level  and  15,000  ft.  A further  comparison 
was  made  with  the  near-f ield/far-f ield  option  (with  low-density  paneling) 
applied  to  the  rigid  model.  (For  a description  of  the  near-f ield/far-f ield 
option,  see  sec.  3.1.3.)  These  correlation  runs  are  listed  as  SD&SS  runs 
1 through  3 in  table  1. 

Having  established  a correlation  of  analyses  between  AFFDL  and  Boeing  for 
the  basic  airplane,  the  balance  of  the  analyses  was  devoted  to  investigating 
modified  paneling  configurations  and  to  the  further  application  of  the  near- 
field/far-f ield  option. 

These  modified  paneling  configurations  are  listed  as  SD&SS  runs  4 through 
8 in  table  1.  Eacli  successive  paneling  configuration  ’as  based  on  the 
experience  obtained  from  the  previous  model,  with  the  objective  being  to 
improve  both  the  hinge  moments  and  airplane  derivatives. 

2.2.1  Comparison  of  AFFDL  and  Boeing  Results  for  Basic  Paneled  Model 

The  basic  paneled  model  analyzed  by  the  AFFDL  is  shown  in  figure  1.  This 
model  was  used  to  compare  directly  the  AFFDL  and  Boeing  results. 

To  check  the  quality  of  the  basic  model,  the  antisymmetric  stability  and 
control  derivatives  for  both  the  rigid  and  flexible  airplane  as  obtained  by 
the  Boelng/AFFDL  FLEXSTAB  analyses  were  compared  (table  2)  with  wind-tunnel 
data  obtained  from  reference  4. 

Table  3 compares  the  hinge  moments  for  the  flaperon,  horizontal  tail,  and 
rudder  obtained  by  the  AFFDL  and  Boeing  analyses  for  the  rigid  airplane  and 
the  flexible  airplane  at  sea  level  and  15,000  ft.  For  the  rigid  airplane  the 
results  (except  for  a few  minor  discrepancies)  are  in  agreement.  However, 
large  differences  are  apparent  in  the  comparison  of  the  hinge-moment  deriva- 
tives for  the  flexible  airplane  both  at  sea  level  and  at  15,000  ft.  The 
rudder  hinge-moment  derivatives  show  the  smallest  differences  between  the 
Boeing  and  AFFDL  results.  The  reasons  for  these  differences  are  not  known, 
but  they  could  be  due  to  using  a different  NASTAP  or  a different  version  of 
ESIC. 
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Table  4 compares  the  AFFDL  and  Boeing  flaperon  hinge-moment  derivatives 
using  the  near-f ield/far-f ield  option.  The  low-density  paneling  case  was  com- 
pared for  the  flaperon  only  because  it  was  evident  from  the  AFFDL  results  that 

while  the  basic  flaperon  hinge-moment  derivatives  C,  , **,  and  ^ compared 
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poorly  with  wind-tunnel  data,  application  of  the  near-field/far-field  option 
considerably  improved  them  as  compared  to  the  wind-tunnel  data. 

While  the  basic  flaperon  hinge-moment  results  between  AFFDL  and  Boeing 
correlated  well,  as  noted  above,  this  is  not  the  case  when  using  the  near- 
field/far-field  option.  The  Boeing  near-field/far-field  results  are  more 
commensurate  with  the  basic  results  (as  is  to  be  expected)  than  are  the  AFFDL 
results  and  generally  show  much  smaller  changes. 

Because  AFFDL  applied  the  neai'-field/far-f ield  option  to  all  the  control 
surfaces  while  Boeing  applied  it  to  the  flaperon  only,  there  are  a different 
number  of  thin  bodies  and  panels  involved.  Also,  it  is  concluded  that  there 
are  some  coding  errors  in  the  option;  this  may  explain  the  differences  in  the 
resuits . 

2.2.2  Investigation  of  Other  Paneling  Configurations 

Having  completed  the  comparison  runs  with  the  basic  AFFDL  paneling,  alter- 
native configurations  were  investigated;  however,  while  this  work  was  in 
progress,  it  was  discovered  that  there  was  an  error  in  the  hinge-line  defini- 
tion of  the  flaperon  in  the  AFFDL  input  data.  Although  the  x-c.oordinates 
describing  the  ninge  line  were  correct,  the  y-coordinates  were  displaced 
inboard,  which  resulted  in  the  flaperon  hinge  line  being  further  aft  than  the 
correct  position.  This  error  is  the  direct  cause  of  the  flaperon  hinge  moments 
being  so  small  compared  to  the  wind-tunnel  data.  As  shown  in  table  5,  cor- 
rection of  the  hinge  line  greatly  improves  the  flaperon  hinge  moments. 

By  the.  same  argument  the  horlzontal-ta  11  hinge  moments  could  be  improved 
since,  compared  with  wind  tunnel  data,  +Cj,  is  too  large  and  -C|,^  is  too 
small.  It  is  evident  that  by  redefining  the  horizontal-tail  hinge  line  a 
little  further  forward,  +Ch  would  be  reduced  and  -C^  increased. 

Table  6 compares  the  hinge-moment  derivatives  obtained  from  the  various 
paneling  configurations  investigated.  In  each  case  the  flaperon  hinge  line 
is  in  the  correct  location.  In  table  6 the  first  model  is  the  basic  paneled 
model  with  the  corrected  flaperon  hinge  line  (fig.  1). 
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The  second  model  has  the  same  paneling  as  the  basic  model  except  that  the 
iAinermost  row  of  panels  along  the  interference  body  has  been  incorporated  into 
one  thin  body  called  STRAKE.  Previously  the  forward  panels  were  part  of  the 
wing  and  the  aft  panels  comprised  a thin  body  called  IHTAIL.  In  general,  it  is 
better  to  have  adjacent  thin  bodies  in  contact  along  chordwise  lines  rather 
than  spanwise  lines,  although  this  is  more  important  when  leading-edge  thrust 
or  thickness  is  being  considered. 

The  third  column  of  tabl~  6 lists  the  hinge-moment  derivatives  for  the 
basic  paneled  model  with  the  low-density  near-f ield/far-f ield  option  applied  to 

the  flaperon  only.  The  flaperon  hinge  moments  ^ and  Cu  ^ and,  especially, 
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are  improved  over  the  corresponding  hinge  moments  of  the  basic  paneled 
a 

model  when  compared  with  the  wind-tunnel  data  in  table  7. 

The  fourth  column  of  table  6 lists  the  hinge-moment  derivatives  for  the 
model  that  was  paneled  to  simulate  the  low-density  near-f ield/far-f ield  model 
(fig.  2).  This  model  gives,  in  general,  no  better  results  than  the  basic 
model  while  being  much  more  expensive  (293  versus  254  singularities). 

The  fifth  column  contains  the  hinge-moment  derivatives  for  an  "improved" 
paneled  model  (fig.  3).  The  paneling  for  this  model  is  changed  from  the  basic 
in  the  following  ways: 

1)  The  wing  panels  are  denser  both  spanwise  and  chordwise  so  that  the 
disparity  in  size  between  the  panels  just  ahead  of  the  flaperon  hinge 
line  and  those  just  behind  is  reduced. 

2)  The  triangular  panels  in  the  basic  model  (beMeved  not  to  be  good 
from  a numerical  point  of  view)  are  removed  and  the  fine  panel!’-  on 
the  flaperon  is  continued  across  the  interference  body. 

3)  The  horizontal-tail  paneling  is  changed  in  order  that  the  st Teamwise 
panel  intersections  continue  to  line  up  with  those  on  the  wing.  The 
paneling  on  the  interference  body  is  lined  up  with  that  on  the 
horizontal  tail. 

4)  The  final  row  of  panels  on  the  interference  body  is  removed  in  order 
partially  to  represent  the  gaps  between  the  empennage  and  the  fuselage. 

5)  The  fin  and  rudder  paneling  is  left  unchanged  since  the  rudder  hinge- 
moment  results  are  good. 


5 


In  spite  of  these  various  modifications  and  a farther  increase  in  the  number 
of  singularities  (299)  over  the  Low-density  simulation  (293),  comparison  to 
wind-tunnel  data  shows  fhrt  the  hinge  moments  are  not  improved  correspondingly. 

SECTION  111 

REVIEW  AND  DISCUSSION  OR  EXISTING  HINGE-MOMENT  ANALYSIS  METHODS 


In  view  of  the  complexity  of  the  subject,  the  discussion  here  is  limited 
to  longitudinal  control  in  the  subsonic  flight  regime.  The  general  discussion, 
however,  can  be  extended  to  lateral -directional  control  and  Lo  the  supersonic 
flight  regime  where  applicable.  The  control  devices  include  various  types  of 
flaps  and  slats,  such  as  plain  or  split  flaps,  single-  or  double-slotted  flaps, 
Fowler  flaps,  leading-edge  flaps,  and  slats.  Complicating  the  analysis  of 
hinge  moments  are  geometric  arrangements  such  as  basic  airfoil  characteristics, 
trailing-edge  bluntness  and  angle,  gap  size  and  geometry,  nose  shape,  and  bal- 
ance ratio.  Except  for  very  small  angles  of  attack  and  control  deflections, 
the  hinge  moments  are  basically  nonlinear.  Al  larger  angles,  the  flow  is 
bound  to  separate  on  the  upper  surface  of  the  trailing-edge  control  surface 
and  hence  is  characterized  by  the  boundary-layer  formulation.  All  of  these^ 
factors  make  the  prediction  of  the  hinge  moment  extremely  difficult. 


A typical  Longitudinal  control  may  be  provided  by  deflecting  Lhe  elevator, 
which  may  be  all  or  part  of  the  tail,  or  a trailing-edge  flap  in  a tailless 
design.  In  addition,  a tab  may  also  be  mounted  on  the  elevator  or  flap.  For 
small  control  deflections,  6,  as  is  usually  stipulated,  one  may  assume  that 
the  hinge  moment,  C|,,  is  linear  in  the  control  deflection  such  that 


C],  - C’i,(u,q)  + C], 
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(1) 

DC), 

x-v— . In  the  case 


where  C^(it,q)  is  the  basic  hinge  moment  when  6 = 0 and  C^  - 
of  small  angle  of  attack,  u,  and  pitch  rate,  q,  the  linear  relation  holds 
among  the  hinge  moment,  the  pitch  rate,  and  the  angle  of  attack  of  the  surface 


that  contains  the  control  surface.  Then, 

ac 

Cl,  - C ’ + u + C,  q + C,  6 
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In  terms  of  the  angle  of  attack,  u,  we  have 
ch  " chQ  + rt  + ch(j  l>  + Ch($  <S 
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where  the  coefficients  are  related  through  tail  downwash,  etc.  (ref.  5),  as 

3C, 
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Hence,  for  the  range  where  linearity  holds,  the  calculation  of  hinge  moments 
reduces  to  the  evaluation  of  ChQ,  Cj-^,  Cj^,  and  separately.  Basically, 

Ci1q  is  the  camber  effect;  and  are,  respectively,  the  hinge  moments  due 
to  unit  angle  of  attack  and  unit  pitch  rate  at  zero  control  deflection;  and 
is  the  hinge  moment  due  to  unit  control  deflection  angle  with  camber, 
angle  of  attack,  and  pitch  rate  set  to  zero.  A typical  load  distribution  on 
the  control  surface  with  and  without  deflection  is  shown  in  figure  4. 


If  the  control  surface  is  not  deflected,  the  evaluation  of  the  coefficients 
» Ch^ * and  is  no  different  from  that  of  lift  and  moment  of  the  basic 
airfoil  for  which  various  techniques,  either  analytical  or  experimental,  are 
available.  The  prediction  of  however,  is  quite  involved.  Physically , 

when  the  control  surface  is  deflected,  the  geometry  of  the  configuration 
becomes  more  complicated,  either  through  the  change  in  the  geometry  of  the 
lifting  surface  or  through  the  creation  of  the  gap  between  the  basic,  airfoil 
and  the  control  surface.  The  actuation  of  the  control  surface  also  induces 
viscous  effects  that  may  alter  the  boundary  layer  ahead  of  the  control-surface 
leading  edge  and,  possibly,  cause  the  flow  to  separate  from  the  control  surface, 
depending  on  a variety  of  geometric  factors.  Mathematically,  potential-flow 
considerations  tell  us  that  the  deflection  of  the  control  surface  includes 
a discontinuity  in  the  planform  and/or  the  slope;  both  factors  introduce 
singular  behavior  along  the  hinge  line,  the  planform  breaks,  and  the  side 
edges.  Aside  from  the  above  complications,  if  the  flow  separation  does  occur, 
as  is  usual  on  the  upper  surface  of  the  trailing-edge  control  with  moderate 
or  higher  angles,  the  existing  analytical  methods  break  down  completely 
since  the  current  theories  for  inviscid  flow  and  boundary-layer  flow  no 
longer  apply.  Consequently,  wind-tunnel  tests  or  semiempir ical  methods  are 
the  last  resort  for  hinge-moment  predictions.  Nevertheless , it  would  be  more 
beneficial  if  a certain  degree  of  confidence  could  be  achieved  through  the 
development  of  capable  analytical  schemes. 

The  object  of  the  current  task  is  to  review  and  discuss  the  theoretical 
and  numerical  characteristics  of  the  FLEXSTAB  method  (ref.  3)  and  a few  other 
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numerical  methods,  with  particular  emphasis  on  the  prediction  of  the  hinge- 
moment  control  derivative,  Cj^. 

'3.1  FLEXSTAB  METHOD 

The  FLEXSTAB  method  for  predicting  stability  characteristics  of  the  air- 
plane is  based  on  the  numerical  scheme  developed  by  Woodward  (refs.  6,  7,  8). 
The  governing  flow  equation  is  the  linearized  potential  flow  with  boundary 
conditions  satisfied  on  the  mean  surfaces;  i.e.,  in  a linearized  sense. 
Numerically,  the  resulting  finite  set  of  integral  equations,  which  correspond 
to  a set  of  boundary  conditions  applied  at  finite  points,  is  reduced  to  a 
system  of  linear  algebraic  equations  with  coefficients  called  the  "aerodynamic 
influence  coefficients"  (A1C).  The  solutions  of  the  linearized  boundary- 
value  problem  are  given  in  terms  of  perturbation  velocities  and  pressure 
coefficients,  from  which  the  forces,  moments,  and  stability  parameters  are 
consequently  derived.  A brief  description  of  the  aerodynamics  in  the  FLEXSTAB 
system  is  given  below,  followed  by  the  hinge-moment  calculations  and  the  dis- 
cussion of  the  numerical  characteristics  of  the  solutions. 


3.1.1  Aerodynamics  in  the  FLEXSTAB  System 

In  the  FLEXSTAB  system,  the  components  in  a general  airplane  configuration 
are  categorized  as  either  thin  bodies  (wing,  tail,  strut,  pylon)  or  slender 
bodies  (fuselage,  nacelle,  external  storage).  Still  another  category  of  imagi- 
nary bodies,  called  interference  bodies,  is  needed  to  serve  as  the  medium  of 
inter  Terence  between  the  thin  bodies  and  the  slender  bodies. 

Since  both  the  flow  equation  and  the  boundary  condition  are  linearized, 
the  principle  of  superpos i tion  applies.  Thus,  for  a body  having  camber  and 
thickness,  the  total  effect  cun  be  considered  as  the  sum  of  the  Individual 
effects  of  camber  alone  and  thickness.  Each  individual  effect  is  described 
below . 


Thin- Body  Lifting  Problem — It  is  known  that  the  lift  is  associated  witli  f Low 
circulation,  which  can  be  realized  on  a sheet  of  vortleltv.  The  lift  of  a 
thin  body  is  simulated  by  placing  horseshoe  vortices  on  the  mean  cumber  sur- 
face. The  distribution  of  the  vortices  is  made  simpler  by  dividing  the  mean 
surface  into  quadrilateral  panels  with  the  vort  icily  strength  on  each  panel 
specified  as  a const  ant.  The  boundary  condition  on  the  solid  surface  is  such 
that  t lie  cumulative  velocity  normal  to  the  surface  is  identical  to  the 
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kinematic  downvash  ip  on  the  boundary  point  (control  point)  of  the  surface. 

V 

Mathematically,  the  vorticity  strength  S on  each  panel  is  determined  from  the 
equation 

ty)  - [AIC]  (SV)  (5) 

where  each  element  of  the  [AIC]  matrix,  the  aerodynamic  influence  coefficient, 
is  precisely  the  normal  component  of  the  perturbation  velocity  induced  by  a 
vortex  panel  of  unit  strength.  Since  the  lifting  pressure  is  directly  related 
to  the  circulation  or  the  vortex  strength,  the  lifting  pressure  at  each  panel 
centroid  location  is  calculated  as 

{C  } - [CPMV]  {SV}  (6) 

V 

where  [CPM  ] is  a diagonal  matrix  with  each  element  equal  to  2.  When  substi- 
tuting equation  (5),  the  lifting  pressure  is  seen  to  be  related  to  the  boundary 
condition  as 

(Cp>  - l LSC ] (7) 

where 

l LSC ] - [CPMV]  [AIC]"1 


Tliin-body  Thickness  Problem — Since  the  combination  of  sources  with  a uniform 
freestreum  yields  stream  surfaces  that  resemble  solid  surfaces,  the  thickness 
effect  of  thin  bodies  is  simulated  by  replacing  the  vortex  singularity  distri- 
bution in  the  lifting  problem  by  source  singularities.  The  thickness  pressure 
on  the  thin  body  eun  be  calculated  without  solving  the  boundary-value  problem 
since  the  source  strength  is  known  to  be  proportional  to  the  surface  slope 
which  lias  been  specified.  In  the  FLEXS'i’Ab  system,  the  source  distribution  on 
eueh  panel  is  either  constant  (CS)  (for  M > 1)  or  linear  (L8)  (for  M < 1). 

In  general,  we  have 


fC 


[CPMCS] 


{scs} 


+ I 


cpmls ] 


(8) 


where  each  element  of  the  [CPM'  ] matrices  is  related  to  the  U- component  of  the 
perturbation  velocity  induced  by  source  singularity  distributions  (constant  or 
linear)  of  unit  strength. 


S.lender-Body  Thickness  Problem — In  a manner  similar  to  the  modeling  of  the 
thin-body  thickness,  the  slender-body  thickness  can  be  simulated  by  distribut- 
ing source  singularity  on  the  mean  centerline.  To  simplify  the  numerical 
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procedures,  a slender  body  is  divided  transversely  into  many  segments.  The 
constant-strength  line  source  is  used  in  subsonic  flow  and  the  linearly  varying 
line  source  is  used  in  supersonic  flow.  At  a boundary  point  of  each  slender- 
body  segment,  che  flow  tangency  condition  is  imposed.  The  source  strength  is 
derived  from  the  tangency  condition,  and  the  pressure  due  to  thickness  effect 
is  computed  as  in  equation  (8).  Since  the  flow  property  is  symmetric  with 
respect  to  the  mean  centerline  axis,  this  problem  is  also  referred  to  as  the 
axial  flow  problem. 

Slender-Body  Lifting  Problem — For  a cambered  slender  body  or  a slender  body 
subject  to  angle  of  attack  or  sideslip,  a lifting  force  or  side  force  can  be 
developed.  Both  phenomena  are  referred  to  as  the  crossflow  problem.  It  is 
known  that  the  derivative  of  the  axial  flow  potential  with  respect  to  the 
radial  coordinate  identically  satisfies  the  crossflow  equation.  The  lift  on 
a slender  body  is  simulated  by  replacing  the  line  source  distribution  in  the 
thickness  problem  by  the  line  doublet  distribution,  as  the  doublet  potential 
is  derived  from  the  source  potential  by  differentiation.  The  doublet  singular- 
ity can  be  quadratically  distributed  (for  M > 1)  or  even  splined  over  three 
centerline  segments  (for  M • 1).  The  strength  of  the  line  doublet  singularity 
on  eoeh  centerline  segment  is  constrained  by  the  boundary  condition  imposed 
at  a certain  point  on  tbe  mean  centerline  of  the  slender  body.  The  resultant 
lifting  pressure  on  each  segment  cun  be  expressed  similarly  to  equation  (7), 
with  the  reference  areas  properly  accounted  for. 

Multiple-Body  Interference  Problem — If  there  are  more  than  one  body,  either  of 
the  same  or  different  category,  in  the  same  flow  field,  interference  effects 
take  place  between  the  bodies.  These  interference  effects  must  be  separately 
accounted  for  and  then  combined  with  the  results  of  the  solutions  due  to 
Isolated  bodies.  The  mutual  interference  among  bodies  can  be  illustrated  by 
a simple  wing-body  configuration.  The  effect  of  the  camber  and  incidence  of 
the  body  on  tbe  wing  is  expressed  in  terms  of  extra  downwash  on  the  mean  wing 
surface  induced  by  doublet  singularity  on  the  b'jdy.  On  the  other  hand,  the 
presence  of  the  camber  and  the  incidence  on  the  wing  also  affect  the  body. 
However,  tills  Influence  is  only  indirect;  i.e.,  it  Is  expressed  as  the  downwash 
on  an  imaginary  shield — the  interference  body — that  wraps  around  L tie  slender 
body.  The  extra  downwash  on  the  interference  body  is  canceled  by  the  normal 
velocity  induced  by  singularities  distributed  on  the  interference-body  panels. 
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These  panels  are  vortex  panels,  as  are  those  on  the  thin  bodies.  The  resultant 
pressure  on  the  interference  body  is  then  transferred  to  the  slender  body  with 
which  it  is  associated.  in  general,  the  pressure  on  the  slender-body  segments 
and  interference-body  and  thin-body  panels  can  be  expressed  as: 

{Cp}  - [APq]  (9) 

where  {^}  includes  incidence  and  camber  on  slender-body  centerline  segments 
and  thin-body  panels,  while  [Apg]  contains  the  ensemble  of  [LSC]  for  each 
body,  with  the  columns  corresponding  to  interference- body  panels  deleted. 

The  contribution  of  thickness  to  the  lift  (buoyancy  force)  can  likewise  be 
included  as 

{CpVnt  - [LSC]  {i|>S}int  (10) 

where  is  the  negative  of  the  flow  incidence  on  interference-body  and 

thin-body  panels  due  to  source  singularities  that  are  employed  to  simulate  the 
volume  of  the  wing  and  the  body,  i.e., 

4S}i,lt:  - [TD]  [SS]  (11) 

with  [TD]  being  the  thickness  downwasli  matrix.  The  resulting  pressure,  which 
includes  the  thickness  effect,  is  then 

{Cp}  - [Ap0]  {«!>}  - {CpVJlnt  (12) 

3.1.2  Computation  of  Hinge  Moment 

The  lifting  pressure  given  in  equation  (12)  contains  the  loading  on  the 
control  surface,  which  is  part  of  the  thin  body.  The  hinge-moment  coefficient 
is  computed  us 


cl>  ■ SCC1;  Ihci  Si  ACI>1 

■ $j%  |dc!  <V  w - (cPV)l"t) 


where,  the  nonzero  elements  of  [dc]  are  the  distances  from  the  centroid  of  the 
control  surface  panels  to  the  hinge  line.  The  elements  in  [Tpp]  are  the  areas 

of  control  surface  panels,  and  S and  C are,  respectively,  the  reference  area 

c c 

and  length  for  the  hinge-moment  coefficient. 

Note  that  the  incidence  vector  { } consists  of  the  effects  due  to  camber 
and  control-surface  deflections  as  well  as  various  motion  variables,  such  as 


angle  of  attack,  sideslip,  pitching,  rolling,  and  yawing  velocities.  Since 
only  a linear  relationship  is  considered,  we  have 

{i/>}  * (<pcl  + Os)  <5  + {i^a}  a + {^3}  3 

+ {ipp}  p + {^q}  q + {t^R}  r (14) 

For  longitudinal  motion,  only  the  incidences  due  to  camber,  angle  of  attack 
a,  pitch  rate  q,  and  the  elevator  type  of  the  control  deflection  5 are  con- 
sidered. With  equation  (14)  in  mind,  the  hinge-moment  derivatives  for  longi- 
tudinal motion  can  be  obtained  from  equation  (13)  as 

Ch0  “ s^T  [dc]  tTFP]  {[V  {^}  ' {GpV)int} 

Gha  " S~C“  ^Ap01j  ^>ot* 

Chq  ‘ Vc7  [dc]  l'rFP]  tApo] 

Cl>6  - [dd  <15) 

3.1.3  Discussion  of  Numerical  Characteristics 

Without  elaborating  on  the  limitations  that  are  imposed  on  the  potential 
flow  solutions,  it  should  be  noted  that  the  linearization  of  the  flow  equation 
and  boundary  condition  takes  on  different  forma  that  vary  from  one  p;  nel 
method  to  another.  in  the  FLEXSTAB  analysis,  the  aerodynamic  representation 
is  such  that  the  boundary  condition  is  also  applied  in  a linearized  sense; 
Ji.e.,  the  boundary  condition  la  satisfied  only  on  the  mean  surface  rather  than 
the  actual  configuration  surface.  Although  the  approximation  of  the  boundary 
condition  is  fully  consistent  with  the  linearized  theory,  it  does  impose  a 
severe  restriction  on  the  control-surface  deflection  as  far  as  the  computation 
of  the  hinge  moment  is  concerned.  Thus,  the  vortex  singularities  used  to 
simulate  the  lift  on  the  control  surface  are  p Laced  on  a mean  surface  as  if 
the  control  surface  were  not  deflected,  and  the  boundary  condition  is  also 
applied  at  the  mean  surface,  with  the  tangency  condition  (local  incidence) 
evaluated  at  the  deflected  configuration.  Furthermore,  the  FLEXSTAB  analysis 
employs  the  velocity  condition,  while  in  the  linearized  theory  it  lias  been 
shown  that  it  is  the  mass  flux,  instead  of  the  velocity  boundary  condition, 
that  satisfies  the  flow  equation.  That  is  to  say,  the  tangency  condition 
implies  that  the  mass  flux  is  tangent  to  the  configuration  surface.  The 
perturbation  mass  flux  differs  from  the  perturbation  velocity  in  the  component 
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lying  in  the  compressibility  axis  direction  in  which  the  mass  flux  component 
is  equal  to  the  velocity  components  multiplied  by  a factor  (1  - M2). 

Although  the  interference  effect  between  all  aerodynamic  components  has 
been  accounted  for  in  the  manner  inherent  in  equation  (9),  the  role  of  the 
thickness  only  manifests  itself  as  the  effect  of  buoyancy.  Consequently,  the 
thickness  of  an  isolated  airfoil,  in  particular  the  trailing-edge  angle  of 
the  control  surface,  has  no  effect  on  either  the  loading  distribution  in 
equation  (12)  or  the  hinge  moment  and  its  derivatives  in  equations  (13)  and 
(15)  since  {Cp^}^nt;  in  equation  (10)  is  identically  zero.  It  is  known,  how- 
ever, that  both  the  control-surface  leading-edge  bluntness  and  the  trailing- 
edge  angle  are  crucial  in  the  hinge-moment  calculation.  The  information  on 
the  thickness  of  the  airfoil  is  lost  through  the  assumption  that,  in  the 
linearized  theory,  the  flow  problem  can  be  separated  into  the  symmetric  prob- 
lem due  to  thickness  and  the  antisymmetric  problem  due  to  angle  of  attack  and 
camber.  Because  of  the  symmetry  of  the  thickness  problem,  the  thickness  of 
the  airfoil  produces  no  net  effect  on  the  loading.  Thus,  a FLEXSTAB  analysis 
concerning  hinge  moments  is  further  limited  to  very  thin  wings  unless  it  is 
corrected  by  the  local  effect  due  to  thickness,  such  as  the  correction  on 
the  local  velocity  induced  by  the  thickness  of  the  airfoil.  Since  the  lift 
is  associated  with  circulation,  and  hence  the  vorticity,  the  constant-vortex 
panel  implies  a constant-pressure  panel.  In  the  linearized  theory,  the 
singular  behavior  of  the  flow  is  realized  in  the  areas  close  to  a subsonic 
leading  edge,  planform  breaks,  control-surface  hinge  lines,  and  side  edges. 

The  resultant  loading  In  these  regions,  therefore,  exhibits  a large  gradient 
in  a certain  direction.  The  effort  to  simulate  this  type  of  loading  with 
constant-pressure  panels  is  rather  stringent.  An  obvious  way  to  achieve 
accuracy  in  the  solution  is  to  increase  the  density  of  the  panels  in  the 
region  of  the  leading  edge  and  planform  breaks,  and  on  control  surfaces.  This 
may  not  be  a good  practice  from  an  economics  viewpoint  and  because  of  the  core 
limitation  of  the  computer.  Thus,  a better  paneling  scheme  must  be  sought  to 
yield  a loading  distribution  with  a certain  accuracy.  The  panel  scheme  should 
be  intimately  related  to  the  planform  geometry  and  the  kinematic  downwash 
distribution.  One  cannot  overemphasize  the  significance  of  the  solution  near 
the  planform  breaks,  hinge  line,  etc.  However,  a poor  algorithm  concerned 
with  planform  breaks  has  been  built  into  the  FLEXSTAB  program.  If  a stream- 
wise  panel  row  spans  either  a leading-  or  trailing-edge  breakpoint  in  the 
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planform,  the  breakpoint  edge  is  automatically  redefined  by  connecting  the  row 
and  planform  edge  intersections  with  a straight  line  (fig.  5).  The  singular 
behavior  of  the  solutio.  ''ear  the  original  breakpoint  is  therefore  smeared  and 
relocated  to  the  neighboring  regions.  The  result  of  this  practice  will  incur 
more  errors. 

An  alternative  method  for  computing  the  aerodynamic  influence  coefficients 

in  the  regions  of  control  surfaces,  wingtips,  wing-body  intersections,  etc., 

has  been  provided  in  the  FLEXSTAB  program.  This  method  is  the  near-f ield/far- 

field  option,  which  has  been  applied  to  analyze  the  hinge  moment  for  the 

YF-16  airplane  in  section  II.  This  method  presumably  yields  more  accurate 

s 

results  in  the  region  where  the  loading  is  changing  rapidly;  however,  it  has 
had  very  little  evaluation  and  is  basically  unproven.  The  mechanism  of  this 
method  is  described  as  follows.  A near-field  region  may  encompass  a few 
panels.  Each  panel  in  the  region  is  equally  subdivided  into  many  near-field 
panels.  The  downwash-singularity  strength  relationship  is  still  described  as 
in  equation  (5),  except  that  the  incidence  on  the  control  point  of  each  sub- 
panel is  set  equal  to  that  on  the  control  point  of  the  far-field  panel,  which 
has  been  subdivided.  Since  each  near-field  region  has  been  isolated  from  the 
global  region,  the  singularity  distribution  and,  hence,  the  loading  on  the 
region  are  identically  the  same  as  produced  by  an  isolated  wing.  The  diagonal 
elements  of  the  aerodynamic  influence  coefficients  based  on  the  near-field 
calculation  are  then  equivalent  to  those  that  would  yield  the  average  loading 
over  the  near-field  region  as  an  isolated  wing.  As  an  illustration,  we 
assign  a single  panel  in  a many-panel  wing  as  a near-field  region  and  then 
find  the  aerodynamic  influence  coefficient  using  the  near-field  method. 

Suppose  we  properly  subdivide  this  panel  into  10  panel  rows  and  5 panel 
columns  and  then  compute  the  near-field  aerodynamic  influence  coefficient 
matrix  of  size  50  bv  50.  The  singularity  strengths  and  hence  the  lifting 
pressure  on  each  subpanel  will  then  resemble  the  loading  distribution  on  the 
panel  as  an  isolated  wing  corresponding  to  the  flow  incidence  ip.  A 
typical  chordwise  load  distribution  of  a lifting  surface  in  subsonic  flow  is 
depicted  in  figure  6(a).  The  method  states  that  the  aerodynamic  influence 
coefficient  of  this  panel  based  on  the  near-field  calculation  is  proportional 
to  the  incidence  divided  by  the  average  loading  of  this  isolated  panel  wing, 

-ip  [~  /AC  • dx]”1  - -i|>/AC 

Ax  p P 
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On  the  other  hand,  the  chordwise  and  spanwlse  downwash  distributions  induced 
by  a constant-pressure  panel  are  shown  in  figures  6(b)  and  6(c),  respectively. 
The  usual  aerodynamic  influence  coefficient  in  the  far-field  sense  is  the  down- 
wash  evaluated  at  the  control  point  of  the  panel.  Whether  the  near-field/ 
far-field  method  can  yield  a more  accurate  result  is  still  an  open  question, 
as  is  the  optimum  location  of  the  control  point  in  the  panel  method. 

Finally,  this  section  concludes  with  numerical  results  generated  by  the 
FLEXSTAB  program  and  comparisons  with  test  data.  The  test  model  is  an  arrow- 
wing wing-body  combination  provided  with  leading-  and  trailing-edge  flaps.  The 
arrangement  of  panels  used  for  the  subsonic  flow  analysis  is  shown  in  figure 
7.  The  chordwise  pressure  distributions  at  three  spanwise  stations  on  the 
wing  with  the  full-span  trailing  edge  deflected  are  shown  in  figure  8.  Also 
included  in  figure  8 are  the  wind-tunnel  data  given  in  reference  9.  By  com- 
parison, the  FLEXSTAB  method  is  shown  to  be  capable  of  predicting  the  overall 
trend  of  the  loading  distribution  qualitatively.  However,  in  this  case  c*- 
least,  the  predictions  are  quantitatively  poor  near  the  hinge  line,  even  at 
the  inboard  sections  of  this  highly  swept,  low-aspect-ratio  wing. 

3.2  DATCOM  METHODS 

The  Datcom  (Data  Compendium)  methods  of  aerodynamic  stability  and  control 
(ref.  10)  as  compiled  by  the  USAF  are  basically  semiempirical.  This  is  the 
only  reference  included  for  the  discussion  of  semiempirical  hinge-moment  cal- 
culations. The  Datcom  methods  for  computing  hinge-moment  derivatives  are 
limited  to  small  angles  so  that  flow  separation  does  not  take  place.  Under 
such  restrictions,  the  hinge  moments  are  expected  to  be  linear.  Additional 
parameters  that  constrain  the  solution  to  the  linear  range  of  the  hinge 
moments  are  the  trailing-edge  angles,  the  airfoil  thickness  ratios  the  nose 
shape,  and  the  balance  ratio.  The  effects  of  seals  and  gaps  are  important 
for  either  plain  or  balanced  controls,  since  the  pressure  distribu: 'On 
changes  according  to  whether  the  gap  is  sealed  or  unsealed.  Additional 
factors  such  as  type  and/or  location  of  the  seal,  the  nose  shape,  balancing, 
and  the  basic  airfoil  shape  must  also  be  considered.  The  Datcom  methods, 
however,  do  not  account  for  unsealed  gaps. 

The  hinge-moment  derivatives  discussed  in  the  Datcom  methods  are  C|, 
and  C',^.  For  either  case,  the  section  derivatives  are  given  first,  based  on 
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the  lifting-line  theory.  The  hinge-moment  derivatives  are  then  computed  from 
section  characteristics,  using  a lifting-surface  correction  to  account  for 
sweep  and  aspect  ratio.  Ti.e  Mach-number  effect  is  accounted  for  by  using  the 
Prandtl-Glauert  rule.  Hence  the  Datcom  methods  are  suitable  only  for  higher 
aspect-ratio  wings  and  lower  Mach  numbers.  They  do  not  provide  the  detail 
required  for  the  load  distributions  on  control  surfaces.  In  this  section,  an 
excerpt  is  given  from  Datcom  concerning  the  calculations  of  hinge-moment 
derivatives.  This  is  followed  by  an  analysis  of  the  hinge-moment  derivatives, 
Cqla  and  C^.,  for  the  YF-16  wing  configuration  using  the  Datcom  methods,  and  a 
discussion  of  the  methods  in  general. 


3.2.1  Method  of  Computing  Hinge-Moment  Derivative  Cj^ 

The  semiempirical  method  of  Datcom  for  evaluating  section  hinge-moment 

derivative  applies  only  to  sealed  trailing-edge  control  surfaces  in  the 

a 

linear  moment  range.  This  derivative  is  computed  as 


C"-  ‘ ‘«V  theory'  ^ '>“ory 

c 

4-2  (C  ) . fl  - — r 


a 


- ') 


(16) 


where  Ci,  /(Ch  ) , and  Ci>  / (C p ) . are  empirical.  The  second  term 

la  a theory  'a  Sx  theory  1 

accounts  for  the  thickness  distribution,  with  <J>  being  the  trailing-edge  angle 
and  T the  thickness  ratio. 


The  nose  shape  and  the  nose-balance  effects  are  given  by  another  empirica] 
relation  as 

(Cu  ),  . 

•‘a  balance  , 


<°h  ) 


a' balance 


Ml 


(17) 
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and  the  compressibility  effect  is  included  via  the  Prand tl-Clauer t rule 

low-speed 


<C'Vx>f 


A - M? 


(18) 


Tills  secLion  hinge-moment  derivative  is 
derivative  for  an  entire  wing  as 
AK  eos  A,,,. 

n _ t.  / A , \ 

'h((  AR  4-  2 cos  A,,,  ' ‘V'Vl 
C/4 


then 


4-  Ac 


used 


h 


<x 


to  compute  the  hinge-moment 
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where  ACu  is  an  approximate  lifting-surface  correction  accounting  for  induced 

a 

camber  effects. 

The  expression  is  directly  applicable  to  control  surfaces  having  constant 
chord  ratios  and  constant  airfoil  sections.  Otherwise,  the  average  value  of 
the  section  characteristics  should  be  used. 


3.2.2  Method  of  Computing  Hinge-Moment  Derivative  Cftg 


a 


The  Datcom  method  of  computing  is  very  similar  to  that  for  Ch  • It 

o ct 

is  also  subject  to  the  limitations  as  discussed  for  the  Ch„  computations. 
The  basic  section  derivative  is  given  as 

Ch6’ 

l ^h§)t]1eory 
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ctV  - I- 


theory 


+ 2 ('C^6)  theory  [ 1 (C £ ) 


-<$  theory 

The  nose  shape  and  the  balance  effects  are  corrected  as 


] (tan  2 ~ T) 


(20) 
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„ i , balance  , 

C^6  1 


(21) 
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and  the  compressibility  is  likewise  included  by  using  the  Prandtl-Glauert  rule. 
The  hinge  moment  is  computed  from  the  section  characteristics  as 

2 cos  Ac/4 

Ch6  = cos  Ac/4  cos  V [(CVm  + U6  (CVm  aYTII 


+ ACh, 


(22) 


where  is  the  two-dimensional  lift-effectiveness  parameter  calculated  either 


as 


6 


(C4), 


(23) 


5 


or,  if  experimental  values  of  the  section-lift  increment  AC^  are  available, 
as 


(AC.,)  , 

, , L exp 

= \ * 
u o 


(24) 
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3.2.3  Numerical  Results  for  YF  16  Wing 

The  computations  of  hinge-moment  derivatives  for  the  YF-16  wing  using  the 
Datcom  methods  are  given  below. 


Wing  Data: 


b = 348.4  in 
S = 40,320  in2 
AR  = 3 

Cc  = 21.85  in 
Sc  = 4,067.6  in2 
Cc/c  = 0.191 


M = 0.9 

qx  = 1,201.43  lb/f t? 

Balance  ratio  = 0 

Tan  A = 0.15893 
HL 

Tan  A cju  = 0.62635 
No  thickness 


From  equation  (6.1.3.1-a)  of  reference  10,  we  have 


V = — ] 

a ^ ha; theory 


(Ch  )f 

a c 


neory 


per  radian 


Next,  compute  the  Reynolds  number, 


p = 0.00237689  slug/ft3 
= 1,005.4475  ft/sec 
c = 131.244  in 

(J  = 3.719x10"^  slug/ft-sec 
PUM  c -I 

R.e  = ^ — = 7.03x10 

From  figure  4.1.1.2-8a  of  reference  10,  we  find 


(C*  ) 


£ a theory 


- 0.96 


Then,  use  figure  6.1.3.1-lla  of  reference  10, 


theory  0,9 


Also  from  figure  6.1.3.1-llb  of  reference  10, 


Thus , 


(Cu  ) - -0.5  per  radian 
na  theory  ^ 


“ -0.45  per  radian 


‘C-.iWirt 


Correcting  for  Mach  number  by  the  Prandtl-Glauert  rule,  we  have 


«V„  ■ -1-332 

Next  apply  above  section  characteristics  to  the  entire  wing  as  follows: 
equation  (6.1.6.1-a)  of  reference  10, 

AR  cos  A ,, 

Cha  “ AR  + 2 cos  + ACha 

- -0.55886  + ACha 


Then,  from  figure  6.1.6.1-19a  of  reference  10,  we  find 


AC,. 


- 0.014 


B2  K2  cos  Ac/2 
Using  figure  6.1.6.1-19c  of  reference  10,  we  have 
b2  - 0.92 

From  section  4. 1.1. 2 of  reference  10,  we  calculate  C^  as 

ct ...  - Hr  ithtt — -i 


6 <C£,,)  theory 

14.523  per  radian 


theory 


Next , 


r'i  “ b72  " °-2405 

% ■ s?!  - o- 775 

From  figure  6.1.6.1-19  of  reference  10,  we  find 
ay 'U  -1.35 

<VnG  “ 2-90 

then,  use  equation  (6.1.6.1-b)  of  reference  10, 

<KaN  (1  - V - <Ka>n0  U - no) 

Ka  ■ — 


0.6975 


Finally, 

ACu  - 0.1106 
na 


From 
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Hence, 


- -0.55886  + 0.1106  - -0.44828  per  radian 
On  the  other  hand,  the  FLEXSTAB  result  for  an  entire  YF-16  configuration  is 
HM^  « -1,375,229  lb-in  per  radian 


or 


HM 


- -0.1855 


Chct  “ q S C 
nl  c c 

The  comparable  wind-tunnel  result  as  shown  in  section  II  (with  Sc  and  Cc 
adjusted)  is 

C,  ( -«  -0.2127  per  radian 

Next,  from  equation  (6.1.3.2-a)  of  reference  10,  the  section  hinge-moment 
derivative  Chg  is  calculated  as 

ch(5 ' 

c 

with  ~ - 0.191  and  C^/ (Cha) theory  “ 0*96,  we  find  from  figure  6.1.3.2-12a 
of  reference  10, 


V 

(<V  theory 


0.97 


whereas  (C^) theory  ®^ven  by  figure  6.1.3.2-12b  of  reference  10  as  -0.92. 


Thus 


Che'  “ -0.8924  per  radian 


At  M = 0.9,  we  have 


<Ch(5>M  * “2.0473 

The  hinge-moment  derivative  for  the  entire  wing  is  given  in  figure 
6.1.6.2-a  of  reference  10  as 

2 cos  Ac M 


C|>6  ■ cos  Ac/4  co“  AHL  [ (<Vm  + “6  (<Vm  STZ 

Now,  the  definition  of  the  lift-effectiveness  parameter  is 


1 ] + ACh. 

cos  Ac/4  h6 


-(C^a  /(%>S 
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R 


where  (C£g)a  can  be  obtained  from  equation  (6.1.1.1-c)  of  reference  10,  for 
small  flap  deflection,  as 

C£<5 

<C*4>“  ' theory'  ^theory 


\ \ 


From  figure  6. 1.1.1- 39a  of  reference  10, 

(CV theory  " 3‘4  per  radian 
and  from  figure  6.1.1.1-39b  of  reference  10, 

C*S 


0.93 


theory 


then 


so  that 


Thus 


(C£(S)a  - 3.162 


a&  - -0.218 


Chfi  " -1.646  + AC|1(S 

Now,  from  figure  6.1.6.2-15u  of  reference  10, 

AC1‘6 

- 0.025 


find 


So  that 


“2 

K<S 

cos  cos 

V. 

3.162, 

U2 

- 0.92,  and 

from 

(K..) 

6 Hi 

m 

1.25 

(K.) 

6 Mo 

m 

2.80 

CKts)n± 

(1 

" Mj.)  " <KcS>  T)o 

(1  - 

no  " ni 


0.5975 


Finally , 


Hence  the  hingo-moment  derivative  is 

C|1(^  - -1.646  4-  0.036372  ■ -1.6096  per  radian 
The  FLEXSTAB  result  fox  the  entire  YF-16  airplane  is 
IlM^  “ -598,490.2  lb-in/radian 
or 

m. 

Gj^  * ■ -0.8075  per  radian 

1 ^1‘  c 'e 

The  comparable  wind-tunnel  result  as  shown  in  section  11  (with  and  C 

adjusted)  is  Cjl(j  “ -0.4592  per  rudian.  It  is  noted  that  botli  the  hinge- 
moment  derivatives  Cha  and  computed  from  the  Datoom  methods  are  about 
twice  as  large  as  the  FLEXSTAB  results.  No  conclusion  could  be  made  since  the 
FLEXSTAB's  prediction  is  based  on  the  untiru  YF-16  configuration.  Also,  the 
Datoom  methods  break  down  for  a low-aspect-ratio  wing  as  on  the  YF-16. 

3.2.4  Discussion  of  Datcom  Methods 

The  Datcom  methods  for  hinge-moment  calculation  are  restricted  to  a linear 
relationship  between  small  angles  and  flow  properties.  The  flow  is  also 
assumed  to  be  attached  and  the  gaps  are  sealed.  The  analysis  is  restricted 
to  plain  trailing-edge  flap  types  of  control  surfaces.  The  hinge-moment 
derivatives  are  computed  from  the  section  characteristics  such  as  the  lift- 
effectiveness  parameter  u,.  as  well  as  C^,  (J^,  and  C^.  Basically,  the 

analysis  employs  the  lifting-line  theory  to  calculate  the  two-dimensional 
hinge-moment  derivatives,  with  an  additional  lifting-surface  correction  to 
account  for  sweep.  Consequently,  the  methods  do  not  account  for  wing-body 
interference  effects  that  may  be  important  for  inboard  controls.  The  analysis 
is  directly  applicable  to  control  surfaces  having  constant  chord  ratios  and 
constant  airfoil  contours  across  the  span;  otherwise,  the  average  section 
characteristics  have  to  be  used.  It  does  not  apply  to  configurations  witli 
wing  cutouts  or  with  control  surfaces  that  do  not  have  ends  parallel  to  the 
plane  of  symmetry.  The  method  does  account  for  such  factors  as  balance-chord 
ratio  and  shape,  trai Ling-edge  angle,  and  bluntness  as  well  as  the  nose  shape. 
Since  the  methods  are  based  on  the  Ilf ting-line  theory,  the  calculations  do 
not  include  the  planform  irregularities  such  as  planform  breaks  where  the 
solutions  exhibit  singular  behavior.  In  summary,  the  semiempirical  methods 
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of  Datcora  are  not  superior  to  the  analytical  methods  so  long  as  they  offer  no 
solution  for  separated  flows  or  unsealed  gaps. 

3.3  RHO  4 METHOD 

RHO  4 is  a computer  prop /am  for  predicting  unsteady  lifting-surface  load- 
ings caused  by  motions  of  leading-  and  trailing-edge  control  surfaces  having 
sealed  gaps  (refs.  11,  12,  and  13).  In  the  regions  of  singularities,  namely 
at  hinge  lines  and  side  edges,  the  discontinuities  in  downwash  ure  removed 
and  handled  separately  by  the  technique  developed  in  references  14  and  15. 

The  residual  downwash,  which  is  fairly  smooth,  is  then  processed  using  the 
mode-kernel  function  approach. 

The  methods  of  solution  for  steady  flow  are  described,  followed  by  a 
discussion  of  the  numerical  results. 

3.3.1  Collocation  Method  of  RUO  4 

The  analyses  of  RIIO  4 are  applicable  to  thin  configurations  having  full- 
span  or  multiple  partial-span  control  surfacus  on  the  leading  or  trailing 
edge  in  subsonic  flow.  The  basic  flow  is  described  by  the  linearized  poten- 
tial-flow equation.  The  numerical  scheme  of  solving  the  flow  equation  can  be 
facilitated  by  uui/  g the  Croen’s  function  method.  The  resulting  integral 
equation  is  of  the  form 

W(x, y)  -ff\ C (C*n)  K(x,y;C,n)  dU'l  (23) 

wliure  K(x,y;^,l|)  is  the  kernel  function,  which  defines  the  downwash  ut  the 
field  point  due  to  the  doublut  singularity  on  the  surfucu  point,  und  AC  (f.,i|) 
is  the  lifting  pressure  related  to  singularity  strength.  The  solution  pro- 
cedures differ  according  to  the  condition  of  the  kinematic  downwash  distribu- 
tion, described  below. 

Solution  Procedures  for  Smooth  Downwash  Distribution1 — The  solution  of  the 
integral  equation  is  obtained  by  assuming  that  the  deflections  of  tine  lifting 
surface  may  be  represented  by  a finite  sum  of  discrete  mode  shapes.  For 
steady  flow, 

* “ En.  (x,y)  q,  (26) 

i 

where  N^(x,y)  is  the  ith  mode  sb  pe  and  q^  is  the  ith  generalized  coordinate. 
Corresponding  to  a particular  mode,  we  have  the  lifting  pressure  Ap^,  then 
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(27) 


Ap  -X)(Ap  ) q 
i 1 1 

Using  the  relation  W - Dz/Dt,  equation  (25)  describes  the  link  between  the  mode- 
shape  slope  and  the  lifting  pressure  of  that  particular  mode. 

The  kernel  function  is  further  divided  into  a nonsingular  part  and  a 
singular  part;  the  latter  contains  the  dipole,  the  inverse  square  root,  and 
the  logarithmic  singularities. 

Upon  integration  in  the  chordwise  direction,  the  singularities  in  the 
integrand  of  the  upunwise  integral  can  be  identified.  In  the  evaluation  of 
the  spanwise  portion  of  the  downwash  integral,  the  integration  interval  is 
divided  into  uubregions  and  appropriate  quadrature  formulas  are  applied  to 
ouch  subregion;  i.e.,  the  square-root  quadrature  formula  applied  at  the  ends 
of  the  interval,  the  logarithmic  quadrature  formula  applied  on  either  side  of 
thu  downwash  station,  and  the  Legendre  quadrature  formula  applied  within  Lhu 
remaining  integration  interval.  Near  the  leading  edge,  an  assumed  pressure 
function  that  is  inversely  proportional  to  the  square  root  of  the  distuncu 
from  tl»e  lauding  edge  is  employed,  which  will  provide  the  least  error  in 
downwash  value  near  thu  leading  edge.  The  above  scheme  is  so  implemented 
that  the  solution  will  not  be  too  sensitive  to  the  collocation  stations. 

The  solution  of  the  integral  equation  is  approximated  by  generating  u 
finite  set  of  downwash  distributions,  the  combination  of  which  satisfies  the 
buundury  condition  at  a finite  set  of  control  points.  Each  downwuoh  distribu- 
tion is  obtained  by  evaluating  the  downwuuh  integral,  using  assumed  polynomial 
pressure  distributions  thuL  conform  to  thu  louding  condition  on  the  plan form 
edges.  Mathematically,  the  louding  function  of  the  jth  mode  has  the  following 
form: 

Ap'(/;.n)  - /(V/2p  ~l/EEamn(1)f(ll°(ii)g(n)(i;,n)  (28) 

m n 

where  f and  g are  certain  trigonometric  Functions  and  umn  are  the  unkown 
coefficients  to  be  determined  from  the  integral  equations. 

Solution  Procedures  for  Discontinuous  Downwash  Distributions — The  above  solu- 
tion procedures  for  predicting  loading  over  a smooth  downwash  field  may  fail 
to  produce  reasonable  results  when  applied  to  configurations  having  discontin- 
uous downwash,  which  may  be  caused,  for  instance,  by  control-surface  defiec- 
Lions.  The  technique  used  to  circumvent  this  difficulty  is  the  introduction 
of  the  residual  downwash. 
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The  residual  downwash  is  the  difference  between  the  kinematic  downwash  and 
the  downwash  produced  by  pressure  distribution  Ap^  obtained  from  the  process 
of  matched  asymptotic  expansion,  which  will  provide  the  proper  discontinuous 
downwash  condition  near  the  hinge  line  and  the  side  edges  of  the  control  sur- 
faces (refs.  14  and  15).  These  pressure  distributions  (inner  solution)  are 
valid  only  in  the  localized  regions.  However,  they  can  be  extended  to  the 
outer  region  since  the  distributions  in  the  outer  region  may  take  any  conven- 
ient form  as  long  us  they  satisfy  the  required  loading  condition  along  the  edge 
of  the  planform.  A unique  solution  will  be  determined  from  the  outer  region 
at  u later  stage.  The  usual  condition  imposed  on  these  loading  functions  is 
that  the  pressure  distribution  diminishes  in  proportion  to  the  square  root  of 
the  distance  from  the  edges. 


The  residual  downwash,  defined  us 

“ro.ldu.1  ' "ki.rn.atic  ' «*•»*«•">  dW"  (29) 

Js  expected  to  ho  fairly  smooth  and  therefore  is  umunduble  to  the  standard 
lifting-aurfuco  solution  process,  i.e., 

(30) 


^residual  “ K(x,y;S,n)  <J£.Ji|]  q^ 

After  solving  for  Apj , thu  resultant  lifting  pressure  becomes 


Ap(r„n)  « pAE(f;,n)  +X) Ap^(fi,n)  q 


j 


.1 


M 


(31) 


Thu  hinge  moment  cun  bu  uasily  calculated  using  the  lifting  pressure  computed 
in  equation  (31). 


3.3.2  Discussion  of  RI!0  4 Method 


The  analysis  in  tills  program  is  based  on  thu  mode-kernel  function  method, 
which  uses  a predetermined  polynomial  type  of  loading  distribution  with  coeffi- 
cients determined  by  satisfying  thu  kinematic  downwash  at  certain  collocation 
points  (standard  flat-plate  boundary  points).  Thu  method  is  applicable  to  thin 
lifting  snrfuces  with  twist  via  local  linearization.  The  thickness  effect  can 
be  included  in  the  downwash  boundary  condition  in  terms  of  local  streamwise 
velocity  variations.  The  control  surfaces  are  assumed  to  have  sealed  gaps, 
and  the  balance  ratios  are  set  to  zero.  The  side  edges  arc  required  to  be 
parallel  to  the  plane  of  symmetry.  The  control  surfaces  symmetrically  located 
on  the  Lifting  surfuces  are  allowed  to  deflect  either  symmetrically  or  anti- 
symmetrically. The  input  mode  shapes  have  to  be  defined. 
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The  special  feature  of  tne  RHO  4 method  is  the  manner  in  which  it  handles 
the  discontinuous  downwaah  as  manifested  in  the  control-surface  deflection. 

The  downwash  discontinuity  is  removed  by  using  the  inner  solution,  which  pro- 
vides the  proper  discontinuous  downwash  near  the  hinge  line  and  side  edges. 

The  residual  downwash  is  therefore  smooth  and  amendable  to  the.  standard  kernel 
procedures.  These  procedures  allow  for  the  accurate  prediction  of  the  singular 
behavior  of  the  flow  properties  near  the  irregular  regions  of  the  planform.  In 
particular,  the  method  is  suitable  for  predicting  the  loading  distribution  on 
the  control  surface  and  hence  the  hinge  moment.  The  method,  however,  may  still 
encounter  numerical  difficulties  In  the  regions  of  planform  breaks.  If  the 
downwash  chord  were  located  at,  or  very  close  to,  the  planform  break  stations, 
the  downwaah  integral  would  exhibit  singular  behavior  duo  to  the  oversimplified 
loading  distributions  used  in  solving  the  downwash  integral.  Whenever  the 
planform  breaks  are  present,  the  apunwise  loading  becomes  discontinuous  at  the 
spanwlse  stations,  which  causes  numerical  problems  and  produces  singularities 
in  the  calculated  downwash  sheet.  Although  it  is  suggested  in  the  method  that 
either  the  downwash  chord  be  placed  at  least  u small  distance  away  from  the 
planform  bruuk  stations  or  the  planform  near  the  breakpoints  be  modified  by 
applying  a small  amount  of  smoothing,  the  ueeuruey  of  the  solution  may  be 
degraded.  Presumably,  one  may  include  curtain  appropriate  loading  functions 
corresponding  to  the  inner  solutions  near  the  planform  breaks  in  the  murmur 
exuetly  parallel  to  that  of  applying  tire  inner  solution  near  the  hinge  lino. 
Since  the  method  deals  only  with  the  lifting  surfaces,  the  wing-body  interfer- 
ence oLfuet  Ui  not  included  in  the.  solution.  The  thickness  effect,  especially 
at  the  nose  and  the  trailing  edge,  can  utiter  the  solution  through  the  change 
of  the  local  velocity,  since  the  boundary  condition  Is  applied  in  a linearized 
sense. 

The  numerical  results  on  the  chordwise  load  distribution  obtained  by  the 
RUG  4 motived  are.  given  in  figures  y and  10  for  a full-span  trailing-edge  con- 
trol surface  at  M * 0.21  and  an  inboard  partial-spun  control  surface  at  M * 
0.60,  respectively,  both  predicted  results  agree  well  with  experimental  data. 
For  a delta-wing  model  shown  in  figure  11,  the  hinge-moment  coefficients  are 
also  computed  and  correlate  favorably  with  measurements  for  a leading-edge 
control  and  a tralling-edge  control.  These  comparisons  indicate  that  the 
pressure  distributions  and  the  hinge  moments  can  be  predicted  with  reasonable 
accuracy  using  the  RHO  4 program,  provided  one  is  aware  of  the  limitations  of 
the  method. 
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3.4 


TEA  230  METHOD 


The  Boeing  aerodynamic  program  TEA  230  (refs.  16  and  17)  was  originally 
developed  to  determine  the  aerodynamic  characteristics  of  fan-in-wing  configu- 
rations using  the  subsonic-potential-flow  theory.  This  program  is  applicable 
to  wings  with  deflected  control  surfaces  and  a wide  variety  of  potential-flow 
boundary-value  problems.  Although  the  solution  is  based  on  the  linearized 
potential-flow  equation,  the  nonlinear  effect  has  been  retained  through  the 
application  of  the  exact  boundary  condition.  Mathematically , the  potential- 
flow  solution  is  obtained  by  reducing  the  problem  to  solving  an  integral  equa- 
tion. Numerical  solution  of  the  integral  equation  is  facilitated  by  the  dis- 
tributions of  constant-strength  source  sheets  over  the  configuration  surface, 
and  distributions  of  vortex  filaments  or  doublets  inside  the  configuration 
surface  and  on  the  wake.  The  source  and  the  doublet  strengths  are  determined 
by  satisfying  a set  of  boundary  conditions  at  finite  points  on  the  boundary 
surfaces.  Once  the  singularity  strengths  ave  known,  the  flow-field  details 
may  be  calculated.  Basically,  the  solution  for  the  subsonic  flow  is  obtained 
in  the  incompressible  flow  domain  with  the  effect  of  compressibility  accounted 
for  by  means  of  Oothert’s  rule.  The  TEA  230  method  is  an  extension  of  the 
vortex-lattice  method  and  provides  for  aerodynamic  representation  of  thick 
wing  configurations.  A boundary-layer  theory  is  also  included  in  the  program 
with  the  purpose  of  indicating  the  boundary-layer  conditions  in  the  areas  that 
are  prone  to  flow  separation. 

The  basic  elements  of  the  potential-flow  theory  as  used  in  the  program  are 
presented  in  detail  because  they  are  also  used  in  the  more  advanced  panel 
methods.  They  are  followed  by  the  description  of  the  numerical  method  and  a 
discussion  of  the  numerical  characteristics  of  the  solutions  pertinent  to  the 
load  prediction  on  the  control  surfaces. 

3.4.1  Potential-Flow  Theory 

A flow  field  that  is  uniform  in  the  ambient  can  be  considered  as  a poten- 
tial flow  except  in  the  boundary  layer  adjacent  to  the  configuration  surfaces, 
the  trailing  vortex  sheets,  the  jets,  etc.  The  potential-flow  equation  c n 
be  derived  from  the  conservations  of  mass,  momentum,  and  energy.  The  result- 
ing equation  is  in  general  highly  nonlinear.  The  total  flow  is  assumed  to  be 
made  up  of  an  onset  flow  and  a perturbation  flow.  The  perturbation  flow,  due 
to  the  presence  of  a disturbance  in  the  onset  flow,  is  considered  to  be  very 
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small  compared  to  the  onset  flow.  Therefore,  the  governing  equation  of  the 
perturbation  flow  can  be  linearized.  By  using  the  Gothert  rule,  the  linear- 
ized flow  equation  can  be  further  reduced  to  the  Laplace  equation.  This  solu- 
tion is  already  known  since  it  also  governs  a wide  range  of  problems  in 
mathematical  physics. 

The  boundary-value  problem  can  be  stated  as  follows.  The  perturbation 
potential  is  governed  by  the  Laplace  equation 

V2<}>  •=  0 (32) 


subject  to  the  boundary  conditions 

1)  Neumann  type  in  which  2p/dn  is  given 

2)  Dirichl.et  type  in  which  p is  given 

3)  Mixed  (Poincare)  type  in  which  a linear  combination  of  p and  dp/dn  is 
given 

Mathematically,  the  solution  to  equation  (32)  can  be  found  by  expressing 
velocity  potential  at  a point  P in  terms  of  (p  and  3<j>/3n  on  the  boundary  of  the 
flow  field  in  the  following  integral  equation  (ref.  18): 

1 rr  1 d<b  1 f f i 3 ,1 


if, 


(33) 


where  3<|>/3n  can  be  identified  as  the  source  strength  and  p as  the  doublet 
strength.  The  area  S (.fig.  12)  need  not  include  the  boundary  at  infinity 
because  the  potential  is  known  to  be  bounded  at  infinity. 

When  the  exterior  of  the  flow  fie.ld  p'  is  considered,  it  can  be  related  to 
the  potential  inside  the  flow  field  (ref.  18)  as 


iis  + h < 7 > dS 


if 


1 341' 

S r 3n  ’ 


dS 


(34) 


with  3/3n'  = -3/3n.  Equations  (33)  and  (34)  combine  to  yield 


1_ 

4 rr 


JY  - ( 

JJS  r v 


, 

3n  3nf 


~ > ds 


/is  <7>  ds 


(35) 


Note  that  there  are  an  infinite  number  of  combinations  of  the  singularity 
strengths  (3cf)/3n  + 3cf>'/3n')  and  (p  - cf>')  that  will  give  the  same  value  of  p at  P 
In  order  to  render  the  system  deterministic,  we  must  either  specify  one  of  the 
two  singularity  strengths  or  provide  a relationship  between  the  two. 
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The  boundary  of  the  potential  flow  can  now  be  identified.  The  outer  flow 
beyond  the  edge  of  the  boundary  layer  is  essentially  irrotational . Moreover, 
since  the  boundary  layer  is  usually  very  thin,  we  extend  the  potential  flow 
field  to  the  configuration  surface,  but  ignore  the  nonslip  boundary  condition 
(the  "outer"  solution).  Since  the  wake  is  essentially  a segment  of  the  flow 
region  filled  with  vortices  shed  from  the  trailing  edge  of  the  wing,  it  should 
be  excluded  from  the  potential  flow  field,  A simplified  treatment  of  the 
trailing  vortices  is  to  imagine  the  wake  as  a thin  vortex  sheet  across  which  the 
normal  velocity  is  continuous.  With  that  in  mind,  equation  (35)  can  be 
written  as 


. i_  //  I ( i£  + Ml  ) ds  + //  (4) 

47 r JJS»  r C 3n  3n'  } dS  4ir  JJSn  U 

D O 

1 // 


* > t: 


( i ) ds 
r 


4tt 


A<J>  f-  ( - ) dS 
w dn  r 


(36) 


w 

where  A<f>  is  the  jump  of  potential  across  the  wake. 

In  order  to  render  the  representation  of  equation  (36)  unique,  we  postu- 
late certain  conditions  on  the  Interior  flow  <J>',  which  is  still  at  our  disposal. 
If  we  put  <j>*  * <J)  and  allow  the  normal  velocity  across  the  surface  to  be  dis- 
continuous, then  equation  (36)  can  be  interpreted  as  a flow  field  set  up  by 
a distribution  of  sources  on  the  surfaces.  On  the  other  hand,  if  we  put 
3<j>'/3n  = 3<f>/ 9n,  and  allow  the  potential  to  be  discontinuous  across  the  surface, 
the  flow  field  is  disturbed  from  the  onset  flow  by  doublet  distributions  on 
the  configuration  surface  and  on  the  wake.  Thus,  the  specification  of  the 
interior  flow  implies  a selection  of  the  type  of  singularity  distribution  to 
be  placed  on  the  surface.  It  should  be  noted  that  unless  a linear  relation- 
ship between  the  two  types  of  the  singularity  strengths  is  specified  (mixed 
boundary  condition) , the  source  and  the  doublet  should  not  be  allowed  to  be 
placed  on  the  same  portion  of  the  boundary  surface. 


3.4.2  Numerical  Method  of  TEA  230 


In  the  present  method,  the  source  singularity  is  used  on  on-'  part  of  the 

surface  S , and  the  doublet  singularity  is  used  on  other  parts  of  the  surface 

S,  and  on  the  wake  S . The  Neumann  type  of  boundary  condition  for  analysis 
d w 

problems  results  in 
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on  the  body  surfaces 


AiJ)  = A(j) 

w 1 1:. 


on  the  wake 


(37) 


The  singularity  strengths  can  be  uniquely  determined  from  the  integral 
equation 
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where 

, U . Oil)'  \ 
o ■»  ( , ) 

On  On 

is  the  source  strength  and 


is  the  doublet  strength. 


In  application,  for  non til ting  configurations  only  the  first  term  on  the 
right-hand  side  of  equation  (38)  needs  to  be  retained.  For  lifting  configura- 
tions, the  second  term  can  be  applied  to  model  the  shear  surface  such  as  fan 
faces,  and  the  third  term  can  be  used  to  simulate  the  trailing  vortex  sheet. 
However,  the  doublet  sheet  starting  at  the  truilitig  edge  will  imply  a con- 
centrated vorticity  at  the  trailing  edge  oi  the  wing.  This  concentrated  trail- 
ing edge  will  euusu  numerical  problems  because  the  source  strength  near  the 
trull  lug  edge  would  have  to  adjust  itself  abnormally  in  order  to  satisfy  the 
boundary  condition  near  the  trailing  edge.  The  above  numerical  difficulty 
can  be  alleviated  by  providing  a doublet  sheet,  on  the  interior  (preferably  on 
the  camber  surface),  as  un  extension  cf  the  trailing  vortex  sheet. 

The  numerical  scheme  to  solve  the  integral  equation,  equation  (38),  is  the 
standard  panel,  method.  The  source  sheet  on  S is  divided  into  quadrilateral 
panels  and  passes  through  the  midpoints  of  the  straight  lines  connecting  the 
corner  points  oi  a quadrilateral,  on  each  of  which  the  source  strength  is 
approximated  as  a constant.  The  doublet  sheet  on  8^  is  similarly  divided  into 
quadrilaterals.  The  doublet  pane  1 need  not  be  planar  because  the  velocity 
field  induced  by  a constant-strength  doublet  panel  turns  out  to  be  the  same  as 
that  of  a ring  vortex  along  its  perimeter  for  which  the  Bfot-Savart  law  can  bo 
used . 
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For  cases  where  there  are  more  boundary  points  than  the  number  of  panels 
(equations),  the  boundary  conditions  must  be  satisfied  in  a least-squares  sense. 
In  this  program,  a boundary  condition  is  applied  at  the  centroid  of  each 
source  panel  and  at  the  average  point  of  each  doublet  panel.  For  each  addi- 
tional condition  to  be  satisfied,  an  extra  control  point  must  be  designated. 

For  example,  the  Kutta  condition  can  be  imposed  at  the  trailing  edge  by 
requiring  that  the  velocity  at  the  trailing  edge  be  directed  in  the  plane 
bisecting  the  trailing  edge.  This  condition  will  control  the  spanwise  varia- 
tion of  the  wake  doublet  strength. 

With  this  type  of  flexibility,  the  doublet  singularity  sheet  of  the 
internal  lifting  system  can  be  placed  in  an  arbitrary  manner.  If  the  near- 
field solution  sucli  as  the  load  distribution  on  the  wing  is  not  of  interest, 
one  may  eliminate  a number  of  boundary  points  on  the  interna]  lifting  system. 
This  is  done  by  arbitrarily  assigning  a certain  pattern  of  vorticity  strength 
(jump  in  potential),  leaving  only  the  Imposition  of  the  Kutta  condition.  The 
lifting-line  theory  implies  a single,  concentrated  vortex  placed  at  the  quarter- 
chord  location.  Although  the  distribution  of  the  bound  vortices  on  the  inter- 
nal lifting  system  has  no  bearing  on  the  far-field  solution,  it  may  impose  a 
burden  on  the  source  distribution.  This  may  result  in  a large  gradient  on 
the  boundary  surface  close  to  the  point  where  a large  isolated  bound  vorticity 
is  placed  on  the  camber  surface.  The  effect  is  further  magnified  for  chin 
wings  where  the  external  surfaces  are  very  close  to  the  camber  surface. 

It  is  ideal  to  distribute  the  bound  vortices  in  a smooth  manner  so  that 
the  source  strength  on  the  external  surface  will  also  be  smooth.  In  figure 
13,  a number  of  concentrated  bound  vortices  are  shown  to  be  located  at  a set 
of  points  on  the  camber  surface.  Each  concentrated  vorticity  corresponds  to 
a jump  in  potential.  The  only  requirement  is  that  the  total  jump  be  equal  to 
the  value  at  the  trailing  edge  (and  hence  the  value  on  the  wake) . 

The  induced  normal  velocity  at  each  singularity  panel  is  evaluated  and  the 
sum  of  all  the  contributions  on  a panel  set  equal  to  the  known  quantity  at  the 
boundary  point  on  the  panel.  This  results  in  a system  of  linear  algebraic 
equations  in  the  matrix  form 

[ AIC ] (A)  - {0}  (39) 

where  f A } contains  |i  and  a. 
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With  source  and  doublet  strengths  solved,  the  flow  properties  such  as 
velocity,  pressure,  forces,  moments,  and  streamlines  can  be  derived  from  equa- 
tion (36).  For  instance,  the  pressure  coefficient  is  calculated  from  the 
Bernoulli  integral  as 


which  is  then  scaled  to  the  value  in  the  compressible  flow  domain  using  the 
Gothert  rule. 

3.4.3  Discussion  of  TEA  230  Method  and  Solutions 

The  panels  have  many  applications  as  basic  building  blocks  in  the  TEA  230 
program.  The  source  panel  can  be  considered  as  impermeable  or  permeable  via 
the  specification  of  the  normal  velocity.  The  permeable  panel  is  used  to 
simulate  the  mass  transfer  through  surfaces,  such  as  controlling  the  inflow  to 
the  fan  face  or  adding  the  displacement  thickness  to  simulate  the  boundary- 
layer  effect.  A boundary  condition  can  also  be  specified  at  an  off-body  point 
for  example,  on  the  engine  inlet.  The  best  result  will  be  obtained  by  having 
a smooth  source  distribution  with  a small  singularity  gradient  on  each  control 
point.  A constant-strength  doublet  panel  is  a basic  element  of  the  lifting- 
system  building  blocks.  Each  doublet  panel  can  be  viewed  as  a ring  vortex. 
Many  constant-strength  doublet  panels  can  be  positioned  end  to  end  to  form  a 
doublet  column.  Doublet  columns  can  be  in  turn  placed  side  by  side  to  form  a 
doublet  lattice  network. 


By  proper  arrangement  of  these  building  blocks,  a polynomial  variation  of 
the  doublet  strength  on  the  network  can  be  achieved.  The  versatility  of  this 
arrangement  can  be  illustrated.  For  instance,  eacli  network  can  be  overlapped 
to  form  a new  element,  thus  providing  an  endless  possibility  for  simulating  a 
variety  of  lifting  systems.  These  networks  can  provide  the  proper  camber  and 
twist  of  a surface  and,  therefore,  can  be  used  in  the  design  of  a surface  for 
a given  load  distribution. 

The  superiority  of  this  method  compared  with  the  FI. EXSTAB  method  or  its 
variants  is  its  ability  to  model  the  actual  configuration  surface  and  the 
applicability  of  the  boundary  condition  in  the  exact  sense.  Consequently,  it 
will  enable  us  to  analyze  the  control-surface  load  distribution  with  provi- 
sions to  account  for  the  nose  shape  and  the  trail ing-edge  angle.  As  a matter 
of  fact,  the  TEA  230  program  is  more  capable  of  analyzing  a thick  wing  than  a 
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cambered  wing  because  of  the  manner  in  which  the  lifting  system  is  implemented 
in  the  program.  For  a very  thin  airfoil,  a high  density  of  internal  panels  is 
needed  to  prevent  strong,  localized  effects  on  the  source  panels.  In  addition, 
since  the  wing-body  interference  effect  takes  place  on  the  actual  configuration 
surfaces,  it  eliminates  the  necessity  for  the  interference  shell  employed  in 
the  FLEXSTAB  system.  This  allows  a greater  flexibility  of  paneling  since  it 
does  not  require  panel  edges  to  be  aligned  in  the  streamwise  direction. 

The  method  of  TEA  230  yields  a pressure  distribution  on  the  source  network 
surface;  i.e.,  on  the  configuration  surface.  In  terms  of  wing  and  tail,  this 
implies  the  pressure  on  the  upper  and  the  lower  surface.  The  load  on  the  con- 
trol surface  can  be  taken  as  the  difference  of  the  pressures  on  the  two  sur- 
faces. Alternatively,  the  hinge  moment  can  be  computed  directly  from  the 
forces  acting  on  each  panel  on  the  upper  and  the  lower  control-surface  source 
network.  For  a load  prediction  near  the  planform  breaks,  the  method  of  calcu- 
lation does  not  differ  from  that  of  FLEXSTAB;  i.e.,  neither  method  allows  for 
the  singular  behavior  near  the  breaks.  For  a control-surface  load  calcula- 
tion, TEA  230  Is  superior  in  that  it  properly  includes  the  thickness 

effect.  The  TEA  230  method  also  applies  the  boundary  condition  on  the 
deflected  control  surface,  although  the  boundary  condition  can  be  applied,  as 
in  the  FLEXSTAB  program,  in  a linearized  sense.  In  essence,  the  linearized 
control  surface  Is  modeled  as  if  It  were  undef looted,  except  that  a local  flow 
velocity  (if  known  a priori)  is  used  to  more  accurately  simulate  the  perturba- 
tion velocity  required  on  the  panel  surface,  i.e., 

v • n “ -U  • n ' V sin  6 (41) 

c x> 

whereas  in  FLEXSTAB,  V is  taken  as  U^.  The  accuracy  of  the  assumption  can  be 
evaluated  by  examining  the  streamlines  resulting  from  the  integration  of  the 
velocity  field. 

The  outboard  wing-pressure  distributions  for  an  aileron  and  an  aileron-tab 
combination  on  a Boeing  747  wing  arc  shown  in  figure  14.  The  solutions  are 
obtained  by  applying  the  boundary  condition  in  a linearized  sense  as  given  in 
equation  (41).  This  approach  was  aimed  at  saving  on  the  computer  costs  since 
only  a reference  configuration  need  be  modeled. 

A comparison  of  the  results  of  this  method  with  experimental  data  is  made 
as  follows  for  a modified  NACA  65-210  airfoil  shown  in  figure  15.  This  wing 
lias  an  inboard  partial-span  tra 1 ling-edge  flap  deflected  at  30  deg.  The 
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computed  spanwise  load  is  shown  in  figure  16  wliiJe  the  forces,  moment,  and 
chordwise  pressure  distributions  at  five  spanwise  stations  are  displayed  along 
with  t he  experimental  data  in  figures  17  and  18.  The  discrepancy  is  most 
likely  the  effect  of  boundary-layer  thickening  or  separation  on  the  flap. 

The  arrow-wing  configuration  (ref.  9)  also  was  analyzed  using  the  TEA  230 
program.  The  panel  arrangement  is  shown  in  figure  19,  and  the  comparison  of 
the  result  with  experiment  is  given  in  figure  20.  Also  included  in  figure  20 
is  the  result  using  the  FLEXSTAB  program.  The  chordwise-lif ting-pressure  dis- 
tribution prediction  by  TEA  230  is  clearly  superior  in  many  respects  to  the 
prediction  of  the  FLEXS'l'AB  method,  especially  in  the  neighborhood  of  the 
trailing-edge  hinge  line.  However,  for  such  a highly  swept,  low-aspect-rat J o 
wing,  botli  potential-flow  programs  are  unable  to  predict  the  flow  characteris- 
tics near  tiie  wingtip  with  a large  flap  deflection.  The  flow  in  that  region 
is  probably  dominated  by  the  spiral  vortex,  which  is  characterized  by  flow 
separation  from  the  leading  edge. 

Since  tiie  TEA  230  method  uses  singularities  of  constant  strungth,  the 
solution  may  also  be  sensitive  to  paneling  of  the  network.  A proper  paneling 
scheme  based  on  loading  conditions  on  tiie  plan  form  breaks  and  discontinuous 
downwash  locations  may  tiave  to  be  sought  so  that  a relatively  sparse  panel" 
arrangement  can  still  predict  the  loading  with  reasonable  accuracy. 

3.5  PAN  Aik  PILOT  PROGRAM  METHOD 

The  FAN  ALK  pilot  code  is  u computer  program  system  that  employs  advanced 
numerical  techniques  in  solving  the  analysis  and  design  problems  tor  steady, 
subsonic  potential  flow  (refs.  19  and  20).  Tiie  outstanding  features  of  tiie 
method  include  curved  panel  approximations,  polynomial  type  of  singularity 
strength  variations,  and  exact  boundary-condition  applications.  The  numerical 
results  are  more  accurate  and  stable  compared  wi tli  other  existing  potential- 
flow  programs.  The  method  is  currently  being  expanded  and  improved  both  in 
engineering  and  software  design  under  a contract  with  several  Government 
agencies  (NASA-ARC,  NASA-LRG , AFFDL) . 

3.5.1  Outline  of  FAN  AIR  Method 

The  Integral  equation  corresponding  to  the  boundary-value  problem  for  the 
potential,  flow  is  the  same  as  that  described  in  section  3.4  for  the  TEA  230 
program.  Specifically,  the  potential  can  be  expressed  as  the  sum  of  the  sur- 
face Integrals 
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For  analysis  problems,  the  spec:  1 f leal  ion  of  the  boundary  condition 
integral  equation 


(42) 


yields  the 
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Auxiliary  conditions  such  as  closure  and  finite  downwash  are  required  lor 
uniqueness  when  solving  equation  (44).  The  sLandurd  panel  schemes  are  used  in 
solving  these  equations.  While  they  are  solved  for  the  incompressible  flow, 
the  subsonic  results  are  subsequently  derived  by  transformation. 

Thu  advanced  technique  for  representing  panels  and  singularity  distribu- 
tions is  designed  to  overcome  certain  practical  difficulties  encountered  in  the 
earlier  panel  methods.  The  flat  panels  witli  constant  singularity  strength  are 
very  crude  in  simulating  boundary  surfaces  having  large  curvature,  sucli  as  the 
wing  leading  edge,  unless  very  dense  paneling  is  used.  The  local  curvature 
and  an  accurate  singularity  strength  gradient  are  also  needed  for  providing  a 
more  realistic,  velocity  gradient  thut  may  be  crucial  in  an  area  such  as 
boundury-luyer  calculations.  Moreover,  most  panel  methods  are  known  to  exhibit 
sensitivity  with  respect  to  the  panel  configuration;  the  combined  features  of 
tills  method  are  to  minimize  this  sort  of  sensitivity.  The  higher  degree  panel 
method  is  required  from  other  considerations  as  well.  For  design  problems,  it 
turns  out  that  a singularity  strength  with  more  than  one  degree  of  freedom  is 
also  needed. 

The  basic  building  block  of  the  paneling  scheme  is  the  network.  A network 
can  be  viewed  us  part  of  the  paneled  configuration  surface  that  is  used  to 
simulate  the  specific  physical  phenomenon  on  that  surface.  These  networks  can 
be  combined  to  ullow  the  analysis  and  design  on  configurations  having  thin  or 
thick  wings,  fuselage,  nacelles,  empennage,  wake,  efflux  tubes,  ducts,  or 
flaps.  Each  network  differs  in  singularity  degrees  of  freedom  and  the  boundary 
points  at  which  the  boundary  conditions  are  enforced.  Four  general  types  of 
the  networks  are  employed  in  tills  method:  source-analysis,  doublet-analysis, 
source-design,  and  doublet-design.  These  networks  are  characterized  by  a 
separate  set  of  the  geometry  specification,  singularity  parameter  location, 
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control  point  location,  and  the  influence  coefficient,  of  course.  The  new 
features  of  the  method  are  briefly  described  as  follows. 


Panel  Surface  Definition — A curve  panel  is  represented  by  fitting  a paraboloid 
through  the  corner  points  that  define  the  configuration  surface.  The  discon- 
tinuities in  surface,  slope,  and  curvature  occur  only  on  the  network  edges. 

Thus,  the  coefficients  in  the  panel  surface  representation 

C ■ Ci52  + can2  + C3£n  + C4C  + Can  + cg  (45) 


are  obtained  by  minimization  of  the  square  error 


R ■“  E \ (r,  - Ck)2  (46) 

k 

where  is  the  ordinate  for  the  ktli  mesh  point  and  is  the  weight  thut  is 
chosen  to  be  very  large  for  a corner  point  through  which  the  surface  will  pass. 
The  ratio  of  the  height  to  the  diameter  | (.  | //(I2  + t)'!  iy  assumed  to  be  small  to 
allow  for  expansion  of  the  kernel  and  the  integrations  in  closed  form. 


Singularity  Strength  Definition — Discrete  values  of  singularity  strength  are 
placed  at  certain  standard  points  on  each  network.  The  distribution  of  singu- 
larity strength  on  a source  network  ia  assumed  to  be  linear,  i.e., 

x(t,n)  - o(e,n)  - aQ  + £ + a n (47) 


and  that  on  a doublet  network  is  quadratic,  i.e., 


x(e,n)  - M(t,n)  - yQ  + K n + \ ^ 

+ 11..  Oi  + \ u 112 
f,n  l mu 


(48) 


The  coefficients  of  the  singularity  distribution  on  each  panel,  source  or 
doublet,  are  obtained  by  minimizing  the  error 

R - yEwk  (A  - Ak)7  (49) 

k 

where  the  summation  extends  over  all  singularity  parameter  locations  on  the 
particular  panel  and  the  neighboring  panels  as  well. 

Boundary  Condition  and  Boundary  Points — At  each  panel,  center  boundary  point, 
the  boundary  condition  for  analysis  problems  reads 
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(50) 


cx(n  • V)  + c2X  - g 
while  for  design  problems 

<4  * V)  + (t2  • VA)  - g (51) 

For  example,  the  usual  condition  for  analyzing  the  flow  over  an  impermeable 
surface  is  given  by  specifying  c. i - 1,  e?  “ g - 0.  In  order  to  control  the  edge 
downwash  or  to  provide  a precise  network  matching,  the  boundary  points  are  ulso 
placed  on  the  network  edges  (which  may  be  slightly  withdrawn  from  the  edges  to 
eliminate  numerical  problems).  Thus,  the  standard  boundary  points  include  the 
panel  center  points  as  rfu.1  1.  as  edge  abutmunt  points  in  case  of  doublet  nutwork 
to  enforce  the  aerodynamic  requirements  such  as  the  Kutta  condition  or  the 
continuity  of  singularity  strength  across  the  network  Junctions,  etc. 

For  the  source-design  network,  it  is  also  required  that  on  uueh  panel 
column  the  closure  condition 


JVle^n  • V)  + c^A]  dS  - // g dS  (52) 

be  satisfied. 

Influence  Coefficient  Calculations — The  potential  due  to  the  singularity 
distribution  on  a curved  panel,  as  it  utunds,  cannot  be  integrated  in  u closed 
form.  However,  the  closed-form  integration  is  superior  to  the  numerical 
integration  in  all  aspects.  The  clouud-form  integrals  are  obtained  for  the 
aerodynamic  influence  coefficients  for  the  near-field  through  an  expansion 
based  on  a small  panel  curvature  assumption;  and  for  the  fur-field,  through 
an  expansion  bused  on  the  large  separation  between  the  singularity  panel  and 
field  point.  Recursion  formulas  that  reduce  the  complicated  expressions  to 
those  appearing  in  the  flat-panel,  constant-strength  singularity  techniques 
are  used  extensively  to  facilitate  the  economy  and  simplicity  in  the  calculation 
of  the  influence  coefficient. 


3.5.2  Discussion  of  FAN  AIR  Method 


In  addition  to  the  various  versatile  aspects  found  in  the  TEA  230  program, 
the  FAN  AIR  pilot  program  lias  many  extra  features. 

With  higher  order  of  singularity  distributions  on  the  panel,  the  accuracy 
of  Lhe  solution  is  presumably  better.  Specifically,  with  linearly  varying 
source  singularity  distribution  and  quadratically  varying  doublet  singularity 
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distribution,  the  computed  perturbation  on  each  panel  is  linearly  varying 
instead  of  constant  as  in  FLEXSTAB  or  TEA  230.  Using  the  Bernoulli  integral, 
the  pressure  is  piecewise  quadratically  varying  on  each  panel.  The  forces  and 
momenta  are  obtained  by  integration  on  each  panel  and  summed  over  all  the 
panels.  The  hinge  moment  calculated  in  this  manner  is  definitely  superior  to 
tlie  constant-strength  singularity  panel  method.  Since  the  spline  techniques 
are  extensively  used  in  the  panel  geometry  definition  and  the  singularity 
strength  distributions,  the  solution  is  more  stable  and  less  sensitive  to 
paneled  configuration.  However,  numerical  difficulties  arise  when  the  udges 
of  some  panels  are  near  the  control  points  of  others,  and  the  use  of  fine 
paneling  in  regions  of  rapid  flow  variation  often  forces  fine  paneling  else- 
whure.  Therefore,  as  in  ail  panel  methods,  the  proper  paneling  criteria  need 
to  be  determined. 

The  numerical  results  are  presented  for  a 35-deg  swoptback  wing  with 
partial-span  t railing-edge,  control  surface  as  depieLed  in  figure  21.  Three 
different  models  are  used  to  simulate  the  flow  on  the  wing,  the  flap,  and  the 
gaps.  In  the  first  two  panel  models,  the  exact  boundary  condition  is  applied 
on  the  flap  surface.  In  the  third  panel  model,  the  ooundary  condition  is 
satisfied  only  in  the  linearised  sense,  similur  to  that  implemented  in  the 
FLEXSTAB  program.  The  modeling  of  the  wake  on  the  side-odge  gaps  differs  in 
the  first  two  panel  models:  for  the  first  panel  model,  it  is  assumed  that  the 
vorticity  is  shed  from  the  trailing  edges  only;  for  the  second  panel  model, 
the  vorticity  is  shed  from  the  w.Lng  side  edges  but  not  from  the  flap  side 
edges.  The  ehordwise  lift  distributions  at  five  spanwise  stations  are 
displayed  in  figure  22.  The  lift  distributions  for  the  first  two  panel  models 
differ  only  along  the  wing  side  of  the  gaps  as  is  expected.  The  difference 
diminishes  toward  the  trailing  edge  and  is  not  discernible  along  the  flap  side 
of  the  gap  (and  lienee  not  shown  in  the  figures).  The  lift  distributions  for 
the  third  panel  model  are  shown  to  agree  better  with  the  FLEXSTAL  results,  as 
they  should.  Both  are  close  to  the  lift  distribution  of  RHQ  4,  which  also 
uses  the  linearized  boundary  conditions.  At  the  midspan,  the  lift  distribu- 
tions corresponding  to  those  three  panel  models  and  the  FLEXSTAB  result  are 
fairly  close.  However,  these  results  indicate  that  the  panel  methods  predict 
higher  lift  than  the  collocation  method  or  the  test  data  given  in  reference  21. 
It  is  unknown  whether  the  discrepancy  is  attributable  to  the  effect  of 
thickness,  which  ia  not  included  in  the  panel-method  analyses. 
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3.6  SUMMARY  OF  METHODOLOGIES  OF  POTENTIAL- FLOW  PROGRAMS 


The  methodologies  of  FLEXSTAB,  RHO  4,  TEA  230,  and  PAN  AIR  pilot  code  ure 
summarized  in  table  8.  The  entries  of  table  8 include  the  year  a particular 
program  was  developed,  the  panel  geometry,  the  singularity  type,  and  the 
restrictions,  if  any.  In  application,  FLEXSTAB  and  RHO  4 programs  only 
accept  mean  surface  paneling;  therefore,  the  deflection  of  control  surfaces  is 
simulated  in  a linearized  sense.  On  the  other  hand,  the  surface  paneling  is 
possible  in  either  the  TEA  230  or  the  PAN  AIR  pilot  program;  lienee,  the  control - 
surface  deflection  can  be  treated  in  an  exact  manner.  However,  associated  with 
the  actual  surface  paneling  is  the  difficulty  of  modeling  the  trailing  vortices 
properly. 


SECTION  IV 

RECOMMENDATIONS  FOR  THE  BEST  APPROACH  USING  THE  FLEXSTAB  METHOD 

Despite  the  numerous  cited  limitations  on  the  FLEXSTAB  method,  the  method 
has  been  found  to  be  udequute  in  predicting  approximate  aerodynamic  character- 
istics, provided:  (1)  the  assumption  of  small  flow  perturbations  is  not 
violated  and  (2)  the  proper  paneling  scheme  is  employed.  For  hinge-moment  cal- 
culations in  the  linear  range,  i.e.t  when  the  control-surface  deflection  is 
small,  condition  (1)  is  likely  to  be  satisfied.  To  satisfy  condition  (2),  how- 
ever, we  need  some  insight  Into  the  locul  aerodynamic  characteristics  along  the 
control-surface  hinge  line,  the  side  edgus,  and  the  plunform  breaks,  inaccura- 
cies arise  in  the  analysis  of  the  loud  distribution  on  a simple  control-surface 
configuration,  not  only  because  condition  (2)  is  difficult  to  achieve  but  also 
because  a constant-pressure  panel  is  inherently  weak  in  handling  the  large 
pressure  gradient  near  the  hinge  line  and  side  edges.  In  this  section,  condi- 
tion (2)  is  explored  in  depth  to  reveal  the  essential  role  it  plays  not  only 
in  the  panel  methods  like  FLEXSTAB  but  in  virtually  any  numerical  method  that 
solves  a differential  equation  using  an  integral  representation.  A better 
scheme  of  handling  the  loadings  in  the  regions  of  singularity  is  also  proposed 
for  panel  methods  in  general. 

4.1  INTEGRAL  REPRESENTATION  OF  DIFFERENTIAL  PROBLEMS 

It  is  known  that  an  equivalence  exists  between  the  integral  equation  and 
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an  initial-value  or  boundary-value  problem  of  mathematical  physics.  In  aero- 
dynamics, these  problems  correspond  to  unsteady  and  steady  flow  problems, 
respectively.  For  fluid  dynamic  problems,  the  partial  differential  equation 
governing  the  conservation  of  mass,  momentum,  and  energy  in  steady  state  is 
linearized  and  then  transformed  into  the  Laplace  equation  for  which  Green's 
function  (or  kernel  function)  is  well  known.  Using  Green's  theorem,  the 
solution  of  the  partial  differential  equation  is  symbolically  expressed  by  an 
integral  equation.  The  solution  of  the  flow  problem  is  then  obtained  by 
solving  the  equivalent  integral  equation  numerically. 

Various  methods  such  us  algebraic,  iterative,  collocation,  and  least 
squares  exist  to  solve  the  integral  equation  numerically.  In  general,  panel 
methods  (FLEXSTAB,  TEA  230,  or  PAN  AIR)  involve  direct  integration  over  the 
panel,  while  the  RHO  4 program  employs  the  collocation  method.  Regardless  of 
the  method  used,  each  results  in  an  algebraic  formulation  of  the  problem 
employing  matrix  algebra,  Two  methods  of  solving  the  integral  equation  are 
reviewed  in  the  next  section  for  identifying  certain  properties  that  may  be 
useful  in  the  panel  methods. 

4.2  METHODS  OF  SOLVING  INTEGRAL  EQUATION 

Consider  the  task  of  determining  A(x)  from  the  linear  integral  equation 

fb  UxiO  A(0  dr  - w(x)  (53) 

*/a 

which  ir  un  integral  equation  of  the  first  kind  or  Fredholm's  equation.  The 
function  K(x;0  represents  a given  kernel  (h*V  operates  on  the  potential 
induced  by  a point  source  or  doublet  singularity),  w(x)  is  a prescribed  func- 
tion (normal  velocity  or  mass  flux  on  the  boundary),  and  A(x)  can  be  inter- 
preted as  the  singularity  strength  or  the  load  distribution. 

With  the  procedures  of  the  algebraic  method,  the  integral  on  the  left  side 
of  equation  (53)  is  replaced  by  a summation  of  n terms  augmented  by  an  error 
term  as 

r\  n 

/ ’ K(x;f)  A (f.)  df  = 2 H.  A(f  ) + E (54) 

J -1  " 

where  the  abscissas  (, . as  well  as  the  corresponding  coefficients  11.  have  been 
chosen  to  obtain  the  desired  degree  of  accuracy  with  a given  number  of  the 
terms  and  to  suppress  t lie  error  term  E below  a predetermined  limit. 
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Inserting  equation  (54)  into  equation  (53)  and  assuming  that  n is  suffi- 
ciently large  and  , 4^  are  properly  chosen  for  the  error  to  be  ignorable, 

the  integral  equation  can  then  be  replaced  by  an  algebraic  relation  of  the 


E Hj  “ w(x)  (55) 

j-1  J 

Since  there  are  n unknown  A(4j)'s  in  equation  (55),  we  can  write  the  same 
equation  for  n discrete  values  of  x as 

E " w(xi>  (56) 

j-1 

i - 1,  2,  . . . , n 

This  represents  a system  of  simultaneous  linear  algebraic  equations  with 
the  same  number  of  unknowns.  The  unknowns  A (4 . ) (j-1, 2, . . . ,n)  can  be  uniquely 
determined  as  long  as  the  determinant  of  the  coefficients  on  the  left  side  of 
equation  (56)  does  not  vanish,  which  is  generally  the  case  for  a well-imposed 
boundary-value  problem.  The  approximate  solution  of  A(x)  in  equation  (53)  can 
be  regarded  as  a polynomial  of  degree  (n-1)  passing  through  the  ri  pivotal 
ordinates. 

In  the  collocation  method,  we  assume  a solution  in  the  form  of  a linear 
combination  of  suitably  chosen  mode  functions  g^(x)  as 

m 

A (x)  ak  gb(x)  (57) 

k-1 

where  a^  are  arbitrary  constants.  Substituting  this  expression  into  equation 
(53),  we  find 


E \f  k(xiO  8k(0  " w(x)  (58) 

k-1  a 

where  the  error  of  representing  A(x)  by  equation  (57)  has  been  ignored. 

When  equation  (58)  is  written  for  ra  distinct  values  of  x^(i-l,2, . . . ,m) , it 
becomes 


E \ MICik  * w<xi> 

k-i  K 
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where 


\ 


MICik  « / K(xi;5)gk(C)  dC 
a 

can  be  evaluated  as  in  equation  (54).  When  the  constants  a (k=l , 2 , . . . ,m)  have 

K. 

been  solved,  the  solution  is  given,  by  equation  (57). 

The  procedures,  as  they  stand,  are  very  simple.  However,  an  error  does 

arise  because  of  ignoring  the  error  term  (apart  from  the  round-off  error).  In 

fact,  the  determination  of  H.  and  which  renders  the  error  term  E a mini- 

J J n 

mum,  encompasses  the  whole  subject  of  the  mechanical  quadratures  (ref.  22). 

Basically,  we  are  concerned  with  a set  of  discrete  abscissas  f, . and  the 

J 

corresponding  coefficients  H^.  , which  will  render  the  equivalence  like  that  in 
equation  (54)  valid  without  the  error  term.  By  Weierstrass'  theorem  on  polyno- 
mial approximation,  any  continuous  function  A(x)  can  be  approximated  within 
[a,b]  to  any  accuracy  by  a polynomial  of  sufficiently  high  degree.  As  a poly- 
nomial of  nth  degree  can  be  uniquely  specified  by  (n+1)  constants,  the  2n  con- 
stants (Hj  and  £ ,j«l,n)  in  equation  (54)  can,  therefore,  accurately  deter- 
mine a polynomial  A(x)  of  degree  (2n-l)  or  less.  If  the  function  A(x)  is  of 
higher  degree,  as  is  usually  the  case  in  aerodynamic  problems,  it  cannot  be 
adequately  represented  by  finite  terms  and  an  error  cannot  be  avoided. 


Graphically,  equation  (54)  without  the  error  term  states  that  the  shaded 
areas  cancel  each  other  between  the  curves;  one  represents  the  (2n-l)st  degree 
polynomial  A(x)  and  the  other  represents  the  polynomial  of  (n-l)st  degree 
that  passes  through  the  nivotal  points,  as  shown  in  figure  23. 

A family  of  the  formulas  that  are  developed  systematically  to  attain  the 
maximum  precision  in  mechanical  quadratures  is  generally  associated  with  the 
Gaussian  quadratures.  Gauss'  work  was  extended  subsequently  by  many  others. 

A detailed  account  of  the  various  types  of  quadratures  can  be  found  in  refer- 
ence 22.  Only  the  essential  elements  of  Gaussian  quadrature  need  to  be 
described  here.  It  has  been  shown  that  the  abscissas  (j=l,n)  chosen  to 
render  the  error  terra  a minimum  (zero  if  A(x)  is  a polynomial  of  (2n-l)st 
degree  or  lower)  are  the  roots  of  a certain  orthogonal  polynomial.  The  type  of 
orthogonal  polynomial  depends  on  the  type  of  the  kernel  and  weight  function  that 
may  be  present.  For  example,  for  a smooth  A(x)  with  unit  weight  function,  the 


called  the  Gauss-Legendre  quadrature;  for  the  logarithmic  type  of  weight 
function  we  use  the  log-Gaussian  quadrature,  etc.  The  properties  and  the 
values  of  the  coefficients  H (j=l,n)  corresponding  to  each  type  of  the  quadra- 

J 

ture  are  also  given  in  reference  22.  The  purpose  of  the  discussion  here  is  to 
shed  some  light  on  the  way  the  panel  methods  should  preferably  be  directed. 

4.3  IMPROVED  FLEXSTAB  METHOD  OF  SOLVING  INTEGRAL  EQUATION 

Without  loss  of  generality,  we  will  consider  the  integral  equation  corre- 
sponding to  the  two-dimensional  problem  given  in  equation  (53).  Assuming  that 

the  function  A(x)  is  expressible  by  its  values  at  n points  as 

n 

A(x)  *=  53  L..(x)  A(x.)  (60) 

3-1  3 

where  L.(j=l,n)  are  obtained  by  certain  interpolation  formulas  or  least-squares 
spline  techniques,  then  equation  (53)  can  be  evaluated  at  x = x^  as  in  equation 
(55),  i.e., 

53  = w<xi')  (61) 

j-1 


where 


H.(x±)  = / Lj(C)  K(Xl;0  d£ 


(62) 


With  the  panel  method,  the  interval  [a,b]  is  subdivided  into  n segments. 
Specifically,  for  the  constant-vortex-panel  method  of  FLEXSTAB,  we  have 

L . (x)  = 1 for  x in  Aj 
J 

- 0 otherwise 

where  Aj  is  the  jth  segment  or  panel.  Then  equation  (62)  becomes 

Hj(xi}  = y*Aj  K(V°  dC  “ AICij  (63) 


so  that  the  solution  of  the  integral  equation  reduces  to  solving  the  matrix 
equation 

[ AI.C  ] { A } = {w}  (64) 

One  notices  that  the  above  equation  is  identical  to  equation  (5)  if  one 
identifies  {A}  as  the  singularity  strength  {s}  and  also  {w}  as  the  downwash 
f'l).  This  linkage  provides  a vehicle  to  convey  the  well-established  knowledge 
developed  in  the  mechanical  quadrature  to  the  panel  scheme  in  the  panel 
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methods.  As  shown  in  section  4.2,  the  values  of  and  f (j«l,n)  must  be  so 
determined  from  K(x;C)  that  equation  (55)  can  be  adequately  used  to  represent 
th<’  integral  equation,  equation  (53),  with  least  error.  However,  the  values  of 
f’.j  are  cxacLly  the  dominant  factors  in  determining  the  paneling  configuration 

for  the  panel  methods. 

The  above  consideration,  however,  Is  incomplete  because  in  the  panel 
methods  thu  control  point  arrangement  is  also  related  to  the  paneling  config- 
uration (even  though  it  need  not  be).  In  general,  the  conLrol  points  at 
(i“l,n),  where  the  boundary  conditions  are  imposed,  are  located  somewhere  within 
the  panel.  Thus,  the  determination  of  •'  will  depend  not  only  on  the  kernel 
function  K(x;0  but  also  on  the  downwasli  distribution  w(x). 


for  instance,  a discontinuous  downwash  field  occurs  whenever  the  control 
surface  Is  deflected,  (hnrespotuling  to  this  downwash  distribution,  there  are 
curtain  loading  character  1st les  along  the  hinge  line  and  side  edges.  The  panel- 
ing configuration  must  therefore  reflect  the  specific,  loading  character  is.  les 
due  to  downwasli  discontinuities.  As  the  downwash  field  is  also  related  Lo  the 
plan  form  geometry,  the  paneling  configuration  is  dependent  upon  the  planform 
shape.  In  general  Llic  loading  characteristics  on  a lifting  surface  depend 
roughly  on  two  distinct  categories  of  configurations.  The  first  category  is 
related  to  planform  edges,  t Ips,  cranks,  and  apexes,  among  which  the  leading- 
edge  m Ingu l./ir i ly  in  subsonic,  flows  and  the  smoothness  of  the  loading  at  trail- 
ing edges  (Kutto  condition)  are  the  familiar  ones.  The  second  category  is 
related  Lo  discontinuous  downwasli  field,  such  as  those  due  to  control-surface 
dell cc i Ions . 


In  reference  23,  the  differential  equation  governing  the  flow  near  the 
planform  breaks  (i.e.,  the  corners,  tips,  cranks,  and  apexes,  fig.  24)  is 
reduced  lo  an  algebraic  eigenvalue  problem.  Kaeh  flow  region  near  the  plan- 
form  break  Is  idealized  us  being  Infinite  in  extent;  i.e.,  an  angular  sector. 
The.  solutions  for  each  region  are  not  unique  because  of  the  lack  ot  a boundary 
condition  at  infiniL’  The  differential  equation  and  the  boundary  conditions 
are  both  homogeneous;  .elutions  can  be  found  only  for  certain  eigenvalues,  and 
the  corresponding  solut  ion  is  an  eigenfunction.  These  solutions  lor  each 
angular  region  arc  the  local  solution’-,  which  are  then  patched  into  the  outer 
solution  valid  on  the  remainder  of  the  flow  field. 


44 


The  Loading  characteristics  near  the  control-surface  hinge  line  and  side 
edges  are  developed  in  reference  14  for  a nonswept  control  surface  using  the 
matched  asymptotic  expansion  method.  The  inner  solution  that  governs  the  flow 
characteristics  near  the  control-surface  leading  edges,  side  edges,  and 
corners  is  further  applied  to  cover  the  special  situation  where  the  control 
hinge  Line  extends  to  the  tip  of  the  lifting  surface,  in  references  11  through 
13  and  15.  A survey  of  the  loading  singularities  inherent  in  the  linearized 
wing  theory  is  given  in  reference  24.  The  available  solutions  are  tabulated 
with  regard  to  the  discontinuities  in  surface  slope,  planform  breaks,  dihedral, 
etc.  The  singularities  are  categorized  as  either  local  or  global.  By  local 
singularity,  the  singularity  strengths — not  only  the  form — are  fixed  by  the 
condition  at  the  singularity  line  and  are  therefore  independent  of  the  general 
shape  and  motion  of  the  lifting  surface.  The  global  singularities  are  those 
that  do  not  permit  the  coefficient  of  singularity  strength  to  be  determined 
prior  to  the  specification  of  the  wing  planform  shape  and  upwash  distribution 
due  to  the  motion.  These  undetermined  coefficients  indicate  the  fact  that  the 
lowest  ord-r  Inner  expansions  of  the  solutions  are  actually  the  eigenvalue  prob 
lem  as  treated  in  reference  23.  In  some  cases,  the  coefficients  can  be  fixed 
through  the  matching  process  where  the  method  of  inner  and  outer  expansion  has 
been  shown  to  apply,  as  illustrated  in  reference  14.  Some  of  the  solutions 
have  been  integrated  into  the  kernel  function  procedures  with  success  and 
Improvement,  as  in  references  11  through  13. 

Once  the  loading  characteristics  near  the  planform  breaks  or  control- 

surface  hinge  line  and  side  edges  are  known,  it  is  possible  to  express  the 

total  loading  I(x)  as  a linear  combination  of  the  singular  loading  functions 

g (x)  corresponding  to  each  particular  region  and  the  regular  loadings  as 
k m n 

*<x>  - £ ak  8k<x)  + £ L-(x>  x'*4>  (65) 

k-1  j=l  J 

Upon  substitutin'.-  into  equation  (53),  we  find  at  x * x^, 
mb  n b 

£ aR  / g (£)  K(x  ;?)  d£  + £ X(C  ) f LAO  KCx^f.)  df.  = w(x  ) (66) 

k=l  a j=l  J a J 1 

Using  the  notation  in  equations  (59)  and  (63),  we  have 
m n 

E \ »mlk  + E 'j  Ate  - (67) 

k-l  J-l 


lifrrrtr  n 
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where  £ w(x^)  and  A(x^). 

Since  there  are  (m+n)  unknowns  on  the  left  side  of  equation  (67)  but  only 
n regular  control  points  for  which  w^(i=l,n)  are  given,  the  system  of  equa- 
tions in  equation  (67)  can  be  made  determinate  by  supplying  m additional  con- 
trol points.  Each  of  these  additional  control  points  is  used  to  determine  the 
coefficient  of  singularity  strength  associated  with  the  singular  region  on  the 
planform  and  is  located  somewhere  in  the  region. 

Recalling  the  methods  of  solving  integral  equations  in  section  4.2,  one 
notices  that  the  solution  procedure  just  outlined  is  equivalent  to  a combina- 
tion of  the  algebraic  method  and  t lie  collocation  method.  The  distinct  feature 
of  the  solution  written  in  the  form  of  equation  (65)  is  that  the  correct  load- 
ing characteristics  near  the  planform  breaks  or  control-surface  hinge  line  and 
side  edges  can  be  properly  accounted  for  by  using  the  accurate  quadrature  inte- 
gration technique  since  the  panel  methods  (especially  ' e constant-pressure 
panel)  are  very  poor  in  handling  the  flow  regions  witli  rapid  change  of  loading 
conditions.  Once  that  part  has  been  taken  care  of,  the  rest  can  be  easily 
handled  by  the  usual  panel  methods. 

In  the  special  case  when  the  singularities  in  these  flow  regions  are  local, 
that  is,  the  singularity  strengths  a^  as  well  as  the  form  g^(x)  are  independ- 
ent of  the  general  shape  and  motion  of  the  lifting  surface  and  they  are  con- 


sidered  as  known  (cf. 

ref . 24) , then 

wc 

can  rewrite  equation  (67)  as 

n 

m 

£ *j 

A1C  . . * w.  - 
ij  i 

E 

a,  MIC., 
k tk 

(68) 

J-l 

k=l 

The  summation  on  the  right-hand  side  of  equation  (68)  is  the  downwash 
associated  with  the  loadings  in  the  singular  regions;  hence,  the  right-hand 
side  becomes  the  residual  downwash,  which  is  presumably  smoother  than  tin* 
original  downwash  distribution.  One  can  find  the  residual  loading  ‘.(j=l,n) 
using  the  standard  technique  in  FLEXSTAB.  The  total  loading  is  the  sum  of  the 
loadings  corresponding  to  the  singular  regions  and  the  residual  loading. 

The  improvement  given  above  is  not  only  applicable  to  the  constant- 
strength  panel  method  FLEXSTAB  but  also  to  other  panel  methods.  In  the  panel 
methods,  the  singularity  strength  on  each  panel  has  been  prescribed  to  he 
either  constant  (e.g.,  FLEXSTAB),  linearly  varying  (e.g.,  ref.  25),  or  of 
higher  order.  In  the  PAN  AIR  pilot  code,  the  doublet  singularity  employed  to 
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simulate  the  lift  is  assumed  to  be  quadratically  varying  and  also  continuous 
across  the  panel  edge  by  a least-squares  spline  fit  over  the  neighboring  panels. 
As  the  degree  of  the  polynomial  assumed  for  the  singularity  distribution 
increases,  so  the  complication  of  the  solution  increases.  As  long  as  the 
degree  of  the  polynomial  for  the  singularity  distribution  is  so  limited,  an 
error  will  exist  that  is  more  pronounced  in  the  singular  regions,  because  the 
loadings  in  those  regions  (e.g.,  logarithmic  function,  inverse  square  root) 
usually  cannot  be  accurately  expressed  by  lower  order  polynomials. 

To  conclude  this  section,  two  numerical  examples  are  given  to  illustrate 
the  effect  of  the  paneling  scheme  on  the  accuracy  of  the  solution.  The  first 
example  is  the  load  distribution  on  the  YF-16  wing  with  flap  deflected  by  0.1 
rad.  The  details  of  the  paneling  of  this  particular  example  are  shown  in  fig- 
ure 1.  The  chordwise  load  distributions  at  three  spanwise  stations  are  shown 
in  figure  25,  (a)  through  (c).  Also  Included  in  the  figure  are  the  residual 
loadings  that  are  obtained  by  subtracting  the  singular  loadings  from  the  total 
loadings  shown  in  the  same  figures.  These  singular  loadings  are  the  inner 
solutions  obtained  from  the  method  of  matched  asymptotic  expansion  (ref.  14), 
with  modifications  applied  to  provide  the  proper  loading  conditions  near  the 
wing  leading  edge,  trailing  edge,  and  side  edges  (refs.  11  through  13).  As 
stated  in  reference  14,  the  inner  solution  of  this  type  yields  not  only  the 
form  but  also  the  strength  of  the  singularity  along  the  hinge  line;  thus  the 
nonvanishing  values  of  the  residual  loadings  near  the  hinge  line  indicate  the 
discrepancy  between  the  predicted  loadings  and  the  loadings  the  linearized 
theory  should  actually  yield.  These  figures  also  show  that  FLEXSTAB  overpre- 
dicts the  loads  over  the  wing  with  control  surface  deflected.  The  downwash 
distribution  corresponding  to  the  singular  loading^  at  each  control  point 
along  three  spanwise  stations  is  shown  together  with  the  step  downwash  distri- 
bution in  figure  25(d).  The  large  value  of  the  residual  loading  along  the 
hinge  line  results  from  improper  paneling  that  does  not  fully  reflect  the 


^ The  downwash  corresponding  to  the  singular  loading  computed  here  is  only 
an  approximation  of  in  equation  (59),  i.e., 

= MICjj  »//K(x1,yi;r,,n)  g^.n)  d£dri 
AE 

' s A,cu  *1  Uj,yi> 

where  (x.,y.)  are  panel  centroid  locations. 
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singular  behavior  near  the  hinge  line  and  the  side  edges.  The  accuracy  of 
the  solution  can  be  improved  if  the  new  schemes,  equations  (65)  through  (68), 
are  used. 

The  second  example  is  a 35-deg  sweptback  wing  of  aspect  ratio  4 as 
depicted  in  figure  21.  To  retain  the  proper  singular  behavior  near  the 
control-surface  hinge  line  and  side  edges  as  predicted  in  the  linearized 
theory,  the  wing  and  partial-span,  trailing-edge  control  surface  are  paneled 
using  the  simple  cosine  function  distribution.  The  total  loadings  and  the 
residual  loadings  along  five  spanwise  stations  are  shown  in  figure  26(a) 
through  (e).  Lt  is  noticed  that  the  residual  loadings  near  the  hinge  line 
(fig.  26(b)  through  (d))  are  much  smaller  than  those  in  the  YF-16  case. 

This  indicates  that  with  proper  paneling  in  the  neighborhood  of  the 
planform  discontinuities,  t ho  FLEXSTAB  results  of  the  airload  on  the  control 
surface  are  reasonably  close  to  those  predicted  by  the  inner  solution  derived 
from  the  method  of  matched  nsymptotic  expansion.  The  downwash  field  corre- 
sponding to  t he  singular  loadings  is  shown  in  figure  26(f)  along  with  a step- 
function  distribution  of  the  prescribed  downwash  due  to  a control  deflection 
of  10  deg.  The  comparisons  of  the  chordwise  pressure  distribution  near  the 
midspan  between  the  FLEXSTAB  result,  the  experimental  data  of  reference  21, 
and  the  collocation  methods  of  references  id  and  26  are  shown  In  figures  27 
and  28.  The  FLEXSTAB  result  is  shown  to  overpredict  the  loading  behind  the 
hinge  line  only  slightly  compared  with  test  data  and  the  results  ol  other 
numerical  methods. 


SECTION  V 

RECOMMENDATIONS  FOR  FUTURE  RESEARCH  AREAS  FOR  H I N(J E-MOMENT  ANALYSIS 

The  subject  of  theoretical  prediction  of  control-surface  hinge  moment 
involves  a wide  spectrum  of  aerodynamic  disciplines.  Although  one  may  con- 
veniently take  the  steady  solution  as  the  lower  limit  of  the  more  general 
unsteady  solution,  the  different  ranges  of  Flight  conditions  (Mach  number,  in 
particular)  make  the  problem  very  complicated  and  thus  beyond  regular  treat- 
ment. For  example,  in  the  transonic  regime,  the  shock  appears  in  the  flow 
field  and  interferes  with  the  boundary  layer.  Under  such  circumstances,  the 
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flow  field  is  nonuniform  in  the  sense  that  the  governing  differential  equation 
in  one  region  has  entirely  different  characteristics  from  the  other.  That  is, 
one  may  have  a supersonic  flow  regime  (flow  equation  being  of  hyperbolic  type) 
embedded  in  the  subsonic  main  flow  (elliptic  type).  In  addition,  the  flow 
equation  in  a layer  adj  :ent  to  the  wing  surface  is  of  parabolic  type.  Tims, 
the  nonuniform  flow  has  deterred  the  general  solution. 

Although  the  major  part  of  the  flow  field  can  be  idealized  as  conservative, 
the  effect  of  viscosity  cannot  be  completely  ignored  (the  Kutta  condition  is 
the  consequence  of  this  effect).  The  effect  of  viscosity  plays  a significant 
role  in  the  upper-surface  flow  of  a positively  deflected  control  surface.  The 

adverse  pressure  gradient  toward  the  trailing  edge  of  the  flap  retards  the  flow 
to  such  an  extent  as  to  cause  the  flow  to  become  separated.  Even  though  claims 

have  been  made  implying  success  in  obtaining  the  approximate  solution  of  the 

flow  phenomena  up  to  the  separation  point  (or  line),  the  solution  beyond  the 
separation  is  still  elusive.  A common  assumption  is  that  the  pressure  is  con- 
stant in  the  bubble  of  separated  and  recirculated  flow.  A reasonable  scheme 
used  to  handle  sucli  a type  of  flow  field  would  be  of  an  iterative  nature.  The 
boundaries  of  different  flow  regimes  (boundary  layer,  thick  wake,  and  potential 
flow)  are  initially  assumed;  then  the  solution  to  each  flow  regime  is  obtained 
until  a solution  is  converged  tliut  satisfies  all  of  the  boundary  conditions. 

The  shape  of  a wake  is  another  aspect  of  the  flow  that  deserves  some  atten- 

tion. Although  the  wake  does  not  signif icuntly  affect  the  flow  on  the  wing 
as  it  does  on  the  tail  or  stabilizer,  a scheme  that  can  estimate  the  self- 
deformation (or  true  position)  of  the  wake  may  prove  to  be  useful  in  imposing 
the  boundary  condition  on  the  aft  lifting  surfaces  and  flaps. 

A new  modeling  technique  for  trailing-edge  wake  analysis  using  a potential- 
flow  program  based  on  the  vortex  method  (Boeing  TEA  372)  was  suggested  in  ref- 
erence 27.  It  was  claimed  that  the  large  fluctuation  in  either  lift  level  or 
lift-curvo*  slope  predicted  by  panel  methods  for  a flapped  wing  arises  from  the 
variation  in  the  path  assigned  to  trailing  vorticity,  especially  that  of  the 
wingtip  and  flap  side  edges.  A reasonably  effective  modeling  scheme  was 
suggested  as  a result  of  helium-bubble-flow  visualization  for  a Boeing 
advanced-technology  airfoil  with  flaps.  Examples  of  modeling  guidelines  are: 
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All  vortices  shed  at  wingtip  should  be  modeled  as  straight  lines  leav- 
ing the  trailing  edge  in  the  freestream  direction. 

2)  The  wake  leaving  the  trailing  edge  of  either  wing  or  flap  should  be 
bent  in  the  direction  of  the  bisector  between  trailing  edge  and  free- 
stream, and  then  aligned  with  the  freestream  after  a short  distance. 

Another  area  of  concern  is  the  flow  phenomenon  near  the  gaps  of  the 
control-surface  side  edges.  It  is  evident  that  the  hingo-moment  calculation 
must  properly  account  for  the  flow  pattern  as  a result  of  the  viscous 
effect.  This  is  the  case  when  mathematical  models  are  devised  to  simulate 
the  actual  flow  problem  by  using  the  compromising  boundary  conditions.  For 
instance,  in  predicting  the  pressure  distribution  on  the  35-deg  sweptback  wing 
using  the  PAN  AIR  pilot  program,  three  different  arbitrary  mathematical  models 
are  used  that  differ  in  the  representation  of  the  gaps  on  the  control-surface 
side  edges  and  in  the  boundary  condition  specifications  (exact  or  approximate). 
The  solution  corresponding  to  each  mathematical  model  differs  significantly  near 
the  side  edges,  as  expected.  It  is  difficult  to  determine  which  model  more 
closely  resembles  the  actual  flow  condition  since  the  supporting  experimental 
data  are  still  lacking.  Once  the  flow  characteristics  are  more  thoroughly 
understood,  an  optimum  model  can  be  constructed  using  a wide  variety  of  the 
options  available  from  the  advanced  panel  aerodynamic  program  currently  being 
developed.  Numerical  experimentation  with  different  mathematical  models  is 
expected  to  be  necessary  in  using  the  potential-flow  program  to  solve  the  real 
flow  problem  as  in  the  deflected  control  surface.  The  concept  guiding  tills 
approach  can  be  Illustrated  by  the  case  of  predicting  load  distribution  on  a 
low-aspect,  highly  swept  wing  set  at  a large  angle  of  attack.  The  load  on  the 
wing  is  dominated  by  separated  flow  right  from  the  leading  edge.  The  potential- 

flow  method  provides  no  means  of  solution  when  the  f Low  becomes  separated. 
However,  from  the  observed  flow  pattern,  the  separated  flow  and  wakes  can  be 
simulated  by  conceiving  a rolling-up  free-vortex  system  with  feeding  sheets, 
fed  sheet,  vortex  core,  etc.  The  boundary  conditions  are  appropriately  applied 
at  the  sheets  whose  locations  are  found  by  iteration  until  the  solution  con- 
verges to  the  one  that  compares  favorably  with  the  measurement  (refs.  28  and  29). 

The  same  concept  can  be  applied  not  only  to  the  flow  modeling  over  the  con- 
trol surface  but  also  to  a wide  variety  of  aerodynamic  problems  such  as  free 
jet,  slipstream,  plume  with  entrainment,  evolution  of  wake,  etc.  In  this 
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sense,  we  can  extend  the  range  of  applicability  of  the  potential  flow  outside 
the  thin  layer  dominated  by  the  viscous  effect  to  include  the  vast  regime  of 
the  flow  field  created  by  the  viscous  effect,  provided  the  mathematical  model 
can  be  devised  to  fit  the  salient  portions  of  the  real  flow  very  well.  Appli- 
cation of  this  sort  is  valuable  since  viscosity  is  always  present.  The  fore- 
going problems  constitute  the  unfinished  aspects  of  aerodynamic  technology  and 
should  be  viewed  as  long-term  research  projects  as  far  as  hinge-moment  predic- 
tion is  concerned. 

The  short-term  research  project  is  improving  the  accuracy  of  the  hinge- 
moment  calculation  based  on  the  linear  aerodynamic  theory.  The  valid  range  of 
linear  theory  is  associated  with  thin,  small-curvature  cambered  wings  at  small- 
\ngle-of-attack  cruising  at  Mach  numbers  away  from  the  critical  value.  Flow 
of  this  type  is  free  of  separation  and  shock  formation,  so  that  equation 
(1)  applies.  Under  such  circumstances,  the  linearized  boundary  condition  may 
be  adequately  applied  since  it  is  consistent  mathematically  with  the  linearized 
flow  equation  and  the  linear  relationship  between  pressure  and  perturbation 
velocity.  All  the  methods  discussed  in  section  III  offer  some  solutions  vary- 
ing to  a certain  extent  from  one  to  another.  The  panel  methods,  in  particular, 
are  superior  because  pressure  loading  functions  need  not  be  assumed  prior  to 
the  solution,  as  required  in  the  collocation  method.  However,  a systematic 
study  of  paneling  configurations  has  yet  to  be  made.  Generally  it  is  not 
possible  to  use  very  dense  paneling  because  of  the  constraints  o ‘ cost  or 
capacity  of  thu  computer. 

In  order  to  alleviate  the  difficulty  of  obtaining  the  optimum  paneling,  which 
is  crucial  to  the  solution  near  the  hinge  line,  side  edges,  or  planforra  breaks, 
the  inner  solutions  obtained  from  the  asymptotic  expansion  method  in  these 
regions  can  be  incorporated  into  the  panel  methods  as  described  in  section  4.3. 
The  recommended  method  subsumes  the  distinctive  features  of  the  existing  numeri- 
cal schemes  in  the  sense  that  the  solution  will  properly  reflect  the  loading 
conditions  in  regions  of  rapid  change  of  flow  properties.  The  philosophical 
aspect  of  the  solution  procedure  is  nothing  more  than  saying  that  the  more  one 
knows  about  the  answer,  the  more  closely  one  can  guess  the  outcome.  The  success 
of  the  proposed  numerical  scheme,  however,  cannot  be  assessed  prior  to  actual 
implementat i on . 
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SECTION  VI 


CONCLUSIONS 

The  success  of  any  method  for  control-surface  hinge-moment  prediction  must 
ultimately  be  evaluated  by  its  ability  to  yield  the  result  that  agrees  with 
experimental  data.  One  criterion  would  be  that  different  approaches  to  a 
mathematical  idealization  of  the  physical  problem  should  yield  the  same  result. 
In  view  of  the  complicated  flow  phenomena  associated  with  the  control  deflec- 
tion, enhanced  physical  understanding  of  the  problem  is  necessary  for  the 
mathematical  treatment  of  the  subject.  The  problem  is  essentially  nonlinear 
in  such  flow  phenomena  as  boundary-layer  separation  evolving  to  free  vortex 
flow  or  shock  waves  interacting  with  boundary  layers.  The  principle  of  super- 
position no  longer  applies  and,  thus,  forms  an  unfinished  area  of  the  subject. 
In  cases  where  the  major  part  of  the  flow  field  can  be  characterized  as  poten- 
tial and  the  perturbation  due  to  the  presence  of  the  body  is  small,  the  linear- 
ized potential-flow  theory  finds  its  application.  Because  of  the  mathematical 
tractability  of  the  theory,  the  elegant  representation  of  the  solution  in  a 
convenient  form  renders  various  numerical  schemes  possible. 

Comparison  of  the  hinge  moments  obtained  by  the  Boeing  and  AFFDL  FLEXSTAB 
analyses  for  the  basic  (AFFDL)  paneled  YF-16  model  showed  the  expected  agree- 
ment for  the  rigid  airplane,  but  the  comparison  for  the  flexible  airplane 
showed  differences  of  up  to  50%  in  some  derivatives.  Since  these  differences 
were  too  frequent  to  be  typographical,  it  is  apparent  that  there  must  be  some 
differences  between  AFFDL  and  Boeing  in  SD&SS  input  data,  NASTAP  data,  or 
the  ESIC  program. 

With  regard  to  the  near-f ield/far-field  option,  the  erratic  results 
obtained  may  be  due  to  coding  errors.  An  investigation  using  a simple  four- 
panel  model  confirms  this  supposition. 

The  importance  of  correct  hinge-line  location  is  obvious  from  the  results 
presented.  In  the  case  of  the  flaperon,  the  hinge  line  was  misplaced  in  the 
AFFDL  results  so  that  the  leading  row  of  panels  subtracted  from,  rather  than 
added  to,  the  hinge  moments,  thus  leading  to  the  very  small  derivatives 
obtained.  A similar  but  smaller  error  may  exist  on  the  horizontal  tail. 
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The  analysis  of  the  four  different  paneling  schemes  employed  (excluding  the 
erratic  near-f ield/ far-f ield  results)  showed  some  variation  between  the  hinge 
moments.  For  the  rudder,  which  is  more  isolated  from  the  other  thin  bodies, 
the  hinge  moments  vary  relatively  little  from  scheme  to  scheme.  However,  for 
the  flaperon  and  (even  more  so)  for  the  horizontal  tail,  the  hinge  moments  are 
more  panel  dependent.  In  the  case  of  the  latter,  this  may  in  part  be  due  to 
the  effect  of  downwash  at  the  tail.  It  is  noticeable  that  the  hinge  moments 
obtained  from  the  low-density  model  were  generally  the  smallest. 

Without  comparing  the  results  with  other  data  (such  as  wind  tunnel  data) , 
it  is  difficult  to  decide  which  paneling  scheme  gives  the  best  results.  Simi- 
larly, there  are  at  present  no  definite  criteria  for  obtaining  the  optimum 
paneling  for  a given  configuration.  Since  the  hinge  moments  are  affected  by 
the  paneling,  it  Is  desirable  that  such  criteria  be  formulated. 

It  is  apparent  from  this  and  other  studies  (e.g.,  ref.  2)  that  generally 
FLEXSTAB  overestimates  the  hinge  moments  C^,  C^,  etc.,  when  compared  with 
wind  tunnel  data.  This  is  because  actual  control-surface  load  distributions 
are  affected  by  such  parameters  as  the  trailing-edge  angle.  FLEXSTAB  does  not 
calculate  such  effects;  however,  the  program  is  useful  in  that  it  calculates 
hinge  moment  and  stability  derivatives  that  are  not  easily  found  in  any  othtr 
way,  together  with  the  ef facts  of  airplane  flexibility  on  their  derivatives. 

Regarding  the  semiempirical  methods  of  Datcom,  the  hinge-moment  prediction 
is  restricted  to  the  linear  range  of  small  angles.  Since  the  method  is  based 
on  the  lifting-line  theory  with  the  effect  of  sweep,  etc.,  accounted  for  by 
using  the  lifting-surface  theory,  its  applicability  is  limited  to  high-aspect- 
ratio  wings  with  sealed  controls.  A method  for  predicting  the  loading  in  the 
transonic  range  is  still  lacking  and  the  wing-body  interaction  is  not  included. 
The  Datcom  method  is,  therefore,  outdated  in  view  of  the  more  sophisticated 
engineering  numerical  methods  available  at  present. 

The  series  or  kernel-function  method  has  the  traditional  beauty  in  that 
only  a minimal  number  of  exact  pressure  functions  are  embedded  in  the  solu- 
tion. However,  a practical  limitation  does  exist  for  selecting  suitable  pres- 
sure functions  to  fit  physical  reality,  such  as  the  loading  at  the  wing-body 
intersection.  Besides  that,  the  truncation  error  is  always  present,  hence 
convergence  appears  to  be  a problem. 
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With  the  advent  of  the  high-speed  computers,  the  number  of  undetermined 
coefficients  need  not  be  so  limited  and  consequently  the  pressure  functions 
need  not  be  assumed.  The  same  problem  can  be  solved  by  the  panel  influence 
methods.  The  panel  method  belongs  also  to  the  finite-element  method  in  t lie 
sense  that  the  boundary  condition  is  only  satisfied  at  finite  points  on  the 
configuration.  With  the  panel  methods,  the  computation  of  the  loading  distri- 
bution on  the  multiple  bodies  (control  surfaces,  in  particular)  poses  no  spe- 
cial problem  from  the  computational  point  of  view,  except  that  the  convergence 
of  the  solution  may  not  always  be  achieved.  In  order  to  reach  a convergent 
solution,  a dense  paneling  near  the  hinge  line  and  side  edges  is  required 
because  the  loadings  in  these  regions  exhibit  singular  behavior  attributable 
to  the  linearized  theory. 

Since  any  computing  machine  is  limited  in  its  capacity  to  handle  the  large 
computational  problem,  a proper  paneling  scheme  that  renders  the  solution  con- 
vergent with  a limited  number  of  panels  needs  to  be  sought.  Thus,  with  proper 
paneling  schemes,  each  numerical  program  may  yield  a solution  that  is  the 
solution  to  the  mathematical  modeling  of  the  physical  problem.  To  enhance 
the  solution  convergence,  a scheme  is  proposed  that  incorporates  the  pressure 
loading  functions  known  from  the  method  of  asymptotic  expansion  near  the 
regions  of  singularity — such  as  those  on  the  planform  leading  edge,  control- 
surface  hinge  line,  and  side  edges — into  the  panel  methods.  The  proposed 
scheme  uses  the  kernel-function  procedures  and  the  panel  methods  so  that  the 
two  categories  serve  to  complement  each  other  to  yield  the  desired  potential- 
flow  solution.  The  value  of  such  a scheme,  however,  can  only  be  measured 
through  further  research  and  application,  which  is  recommended. 
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Figure  3.  YF-16  Mode!  With  "Improved"  Paneling 


Figure  5.  Breakpoint  Edge  Redefinition 


(c)  Spanwise  downwash  distribution  of  a 
constant-pressure  panel 


(b)  Chordwise  downwash  distribution 
of  a constant-pressure  panel 


Figure  6.  Typical  Load  and  Downwash  Distributions  on  a Constant-Pressure  Panel  and  on  a 
Wing  in  Subsonic  Flow 


Figure  8.  Load  Distributions  on  the  Arrow-Wing-Body  Configuration— FLEXSTAB  Versus  Experiment 


Figure  9,  Theoretical  and  Experimental  Chordwise  Load  Distributions  for  a 
Deflected  by  6 * JO  deg  at  a =0  deg  and  M - 0.21 


Full-Span  Flap 
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Figure  10.  Theoretical  and  Experimental  Load  Distributions  for  Partial-Span  Flap 
Deflected  by  b = 10  deg  at  a •*  0 deg  and  M = 0.60 
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pressure  station 


Figure  14.  Outboard  Boeing  747  Wing  Pressure  Distributions  for 


Figure  15.  Source-Pane!  and  Internal  Vortex  Arrangement  on  NASA  Model  Wing  With  Deflected  Flap 
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Figure  }/.  Comparison  of  Theoretical  Force  and  Moment  Hew  Its  With  tx/wrimental 
Data  on  NASA  Model  Wing 
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Figaro  26.  Load  and  Do  wi  wash  Distributions  on  a 36-don  Swap!  back  Win//  With  Control  Surf  tun 
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Figure  27.  Theoretical  and  Experimental  Pressure  Distributions  for  Partial-Span 
Flap  Deflected  by  S m 10  deg  at  a - 0 deg  and  M « 0.60 
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O aw  "0,  6 = 10  deg  (ref.  21, 1?  » 0.46)  — 

□ -0,  S - -10  deg  (ref.  21,rj -0.46)  — 

RHO  4 (tj- 0.4457  ) 0.197 


Figure  28.  Comparison  of  Experimental  and  Calculated  Liftiruj-i'rossuro 

Distributions  on  a Sweptwing  With  Partial-Span  Flap  at  M ~ 0.60 
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TABLE  1 

SD&SS  RUNS  FOR  CORRFLA  T/ON  AND  HINGE-MOMENT  IMPRO  VEMENT  ( ITEM  ZO) 


Run  No. 
1 


Paneling  scheme 

Eiastk 

option 

Attitude 

Near/far 

option 

Basic 

Residual  •elastic 

SL 

No 

Basic 

Q«*kkj  >1 -elastic 

15,000  ft 

No 

Bask 

Rigid 

tb,000  ft 

On  flaperon 
only 

Bask 

Rigid 

SL 

No 

Bask  and  strake 

Rigid 

SL 

No 

Bask 

Rigid 

SL 

On  flaperon 
only 

Low  density 

Residuai-elastk 

SL 

No 

Optimum 

Rigid 

SL 

No 

Notes 


Correlation  with  AFFDL 

Flaperon  hlnge-lln# 
location  incorrect  to  comparr 
with  AFFDL  results 


Redefine  certain  thin  bodies 
to  improve  model 


To  simulate  near-fiald/ 
far-fieid  model 


hinge  location 


TABLE  2 

ANTISYMMETRIC  STABILITY  AND  CONTROL  DERIVATIVES  COMPARISON  OF  FLEXSTAB 
AND  WIND  TUNNEL  DATA  (FROM  REF.  4},  BASIC  PANELED  MODEL 


Wind  tunnel 


Derivative 
(par  red) 


Rigid 


-1.243 


- .1157 


.171 

.014 

- .004 

- .325 
0 

- .038 

- .1339 


[ F textile 

SL 

-1.112 

-1.121 

.239 

.256 

- .086 

.78 

.83 

- .40 

- .428 

.135 

.165 

.013 

.014 

- .003 

- .26 

- .28 

0 

0 

- .036 

- 0356 

- .0770 

- .0922 

.1050 

.126 

- .0673 

- .079 

.0180 

.022 

FLEXSTAB  (Boeing) 

Flexible 

Rigid 

SL 

-00221 

-0.7363 

-0.7637 

.1919 

.1462 

.1619 

- .0761 

- 0796 

.9762 

.8352 

0813 

- .4667 

- .4609 

- .4181 

.1040 

.1171 

.1301 

- .0295 

. .0119 

- .0111 

.0213 

.0173 

.0191 

- .3423 

- .3924 

- .4013 

.1637 

.0994 

.1203 

- .0637 

-.0361 

- .0422 

- .2188 

- .1494 

- .1770 

.4115 

.2441 

.299 

- .2249 

- .1360 

- .1861 

.0679 

.0449 

0476 

10269  deg 

1.3269  deg 

1.9231  dog 

TABLE  3 

COMPARISON  OF  HINGE-MOMENT  DERIVATIVES 
PANELED  MODEL 


-AFFDL  AND  BOEING,  BASIC 


laperon 


AFFDL 

Boeing 

Flexible 

Flex 

ibie 

(per  rad) 

Rigid 

15,000  ft 

SL 

Rigid 

15,000  ft 

SL 

Cha 

-0.000367 

— 

-0.000367 

— 

C*q 

- .00230 

— 

- .00230 

— 

\S 

- .000434 

— 

— 

- .000434 

— 

— 

- .000045 

— 

— 

- .000045 

— 

— 

Ch« 

- .52434 

— 

— 

- .52424 

— 

— 

ChQ 

- .004019 

— 

— 

- .004019 

— 

ch6s 

- .000015 

— 

— 

- .0000147 

— 

— 

ctyj 

.000055 

0.000064, 

0.000071 

.000055 

0.000046 

0.000041 

- .00009 

- .000067 

- .000048 

- .00009 

- .000141 

- .000182 

- .000027 

- .000042 

- .000053 

- .000027 

- .000016 

- .000016 

V 

- .000249 

- .000242 

- .000237 

- .000249 

- .000141 

- .000076 

- .004092 

- .003921 

- .004011 

- .004092 

- .003988 

- .004177 

CK 

- .000324 

- .000312 

- .000305 

- .000324 

- .000302 

- .000287 

Chf 

.000063 

.000061 

.000060 

.000063 

.000061 

.000060 

•V 

- .000013 

- .000012 

- .000012 

- .000013 

- .000012 

- .000010 

t,  I 

■:  t 

H TABLE  3 (CONTD) 


COMPARISON  OF  HINGE-MOMENT DERIVATIVES-AFFDL  AND  BOEING, 
BASIC  PANELED  MODEL 


AFFDL 

Boeing 

Control 

turface 

Derivative 
(per  rad} 

Flexible 

Flexible 

Rigid 

ESSE! 

SL 

Rigid 

1^331 

SL 

Horizontal 

Qic 

0.001921 

0.001921 

_ 

tail 

** 

- .02256 

- 

- 

- .02255 

- 

- 

<*ss 

- .004901 

- 

- 

- .004901 

- 

- 

Chfi 

- .000848 

- 

- 

- .000891 

- 

-8.06998 

- 

- 

-8.05993 

- 

- 

- .106008 

- 

- 

- .106006 

- 

- 

ChfiS 

- .000144 

- 

- 

- .000144 

- 

- 

.001621 

0.002049 

0.002433 

.001621 

0.002283 

0.002857 

.000605 

.000909 

.001357 

.000604 

.001071 

.001547 

- .C ‘SS869 

- .G030  ( . 

- .003258 

- .002643 

- .003145 

- .003536 

Ch  A 

- .003840 

- .004670 

- .005444 

- .003846 

- .005278 

- .006518 

s 

.223056 

.195740 

.181482 

.223054 

.176564 

.144711 

4 

- .004839 

- .005186 

- .005513 

- .004838 

* .005299 

- .005722 

Ch? 

• < % « • ;■ 

• .003367 

- .003897 

- .003001 

- .002242 

< 

- .000127 

- .0001192 

- .000106 

- .000127 

- .0001113 

- .000090 

Rudder 


AFFDL 


Boeing 


Control  Derivative 

surface  (per  rad)  Rigid 


Flexible 


.000132  .000099 


.000412  .000398 


Flexible 

15,000  ft 

SL 

0.000774 

0.000826 

.000094 

.000062 

- .001559 

• .001558 

- .001855 

- .001731 

.066945 

.06235 

.000392 

.000382 

- .001262 

- .001102 

- .000008 

- .000007 

TABLE  4 


COMPARISON  OF  BOEING  AND  AFFDL  FLAPERON  HINGE-MOMENT 
DERIVATIVES  USING  THE  NEAR-FIELD /FAR-FIELD  OPTION 
(LOW-DENSITY  PANELING)  ON  THE  BASIC  MODEL,  RIGID  AIRPLANE 


Deriviative 
(per  rad) 

AFFDL 

Boeing 

Basic 

Near/far 
on  basic 

Basic 

Near/far 
on  basic 

Qna 

• .000367 

- .000867 

- .000367 

- .000483 

- .00230 

- .002603 

- .00230 

- .002235 

chfiS 

- .000434 

- .003812 

- .000434 

- .0004802 

<VU 

- .000045 

- .000136 

- .000045 

- .00015 

- .52434 

.044515 

- .52424 

- .317428 

Ch^ 

- .004019 

- .003542 

- .004019 

- .002636 

- .000015 

- .000005 

- .0000147 

- .0000083 

Ch(j 

.000055 

.000062 

.000055 

.000045 

- .00009 

- .000209 

- .00009 

- .000124 

Chf 

- .000027 

- .000001 

- .000027 

- .000032 

<V 

- .000249 

- .003513 

- .000249 

- .000351 

- .004092 

- .004757 

- .004092 

- .003753 

Chp* 

- .000324 

- .000156 

- .000324 

- .000264 

Ch? 

.000063 

.000074 

.000063 

.000056 

- .000013 

- .000005 

- . W001 3 

- .000009 
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TABLE  5 

COMPARISON  OF  FL  APE  RON  HINGE-MOMENT  DERI  VA  Tl  VES  OF  BASIC  PANEL  ED 
MODEL  WITH  INCORRECT  AND  CORRECT  HINGE  LINES,  RIGID  AIRPLANE 


Control 

surface 

Derivative 
(per  rad) 

Hinge  line 
incorrect 

Hinge  line 
correct 

Flaperon 

cha 

-0,000367 

-0.001146 

Chq 
_ S 

- .002300 

- .006373 

ch<5 

- ,000434 

- .005046 

Chu 

- .000038 

- .000050 

C^a 

- .553539 

-1.48630 

Ch^ 

- .004019 

- .009279 

_ -S 

Chg 

- .0000140 

- .0000306 

ch(j 

.000055 

.000122 

Ch£ 

- .00009 

- .000307 

r< 

- .000030 

- .000080 

Ch6A 

- .000249 

- .004548 

Ch  p 

- .004321 

- .010653 

C^P 

- .000324 

- .000993 

Chf 

A 

.000063 

.000153 

^ 6 

- ,000012 

- .000027 

91 


TABLE  6 

COM  PA  BISON  OF  HINGE-MOMENT  DEB  I VA  77  VES 
FOB  VABIOUS PANELING  CONFIGUBAT/ONS 


Basic  paneling 

^ l-D  NF/FF 

with  strake  on  flaperon 


LD 

paneling 


-0.00106 


Improved 

paneling 


-0.001065 
- .006351 


3 

I 


TABLE  6 (CONTD) 

CO  MPA  RISON  OF  HINGE-MOMENT  DERI  VA  Tl  VES  FOR  VA  RIOUS  PANEL  ING 
CONFIGURATIONS  


Control 

surface 

Derivative 
(per  rad) 

Basic 

paneling 

Basic  paneling 
with  strake 

Basic  paneling 
LD  NF/FF 
on  flaperon 

LD 

paneling 

Improved 

paneling 

Horizontal 

Cha 

0.001921 

0.001933 

0.001493 

0.002485 

0.001755 

tail 

Ch$ 

- .02255 

- .022574 

- .021331 

- .017558 

- .018783 

Ql£S 

- .004901 

- .004895 

- .004888 

- .003062 

- .004068 

Chfl 

- .000998 

- .000939 

- .001014 

- .000891 

- .000901 

04 

-8.51040 

-8.594526 

-7.34355 

-7.882348 

-7.107366 

Ch$ 

- .106008 

- .106329 

- .106233 

- .097523 

- .094317 

Ch6S 

- .000136 

- .000136 

- .000136 

- .000138 

- .000120 

.001621 

.001616 

i 

.001581 

.000965 

.001262 

.000605 

.000608 

.000481 

.000641 

.000396 

Chf 

- .00261 

- .002624 

i 

- .002573 

- .001879 

- .002113 

Ch  A 

- .003846 

- .003844 

- .003844 

- .0022214 

- .0033692 

Ch| 

.235522 

.235626 

.240377 

.249017 

.212466 

< 

- .004839 

- .004858 

- .004516 

- .00444 

- .004311 

Chf 

- .003898 

- .003900 

- .003958 

- .004112 

- .003557 

Ch5*A 

- .0001219 

- .0001219 

- .0001219 

- .0001219 

- .000106 

TABLE  6 (CONCLUDED) 


COM  PA  RISON  OF  HINGE-MOMEN  T DERIVATIVES  FOR  V ARID  US  PA  NEL ING 
CONFIGURATIONS 


Control 

surface 

Derivative 
(per  rad) 

Basic 

paneling 

! 

Basic  paneling 
with  strake 

Basic  paneling 
LD  NF/F 
on  flaperon 

LO 

paneling 

Improved 

paneling 

Rudder 

Ch/3 

0.000700 

i 

0.000700 

0.000701 

0.000637 

0.000651 

% 

.000132 

0.000129 

0.000133 

0.000121 

0.0001 33 

Ch? 

- .001573 

-0.001573 

-0.001573 

-0.001443 

-0.001487 

ch6 

- .002043 

-0.002043 

-0.002043 

-0.001797 

-0.001983 

Ch? 

.083185 

0.083172 

0.083136 

0.075232 

0.080189 

Ch£ 

,000412 

0.000415 

0.000407 

0.000377 

0.000407 

Chf 

••  .001504 

-0.001504 

-0.001503 

-0.001337 

-0.001459 

C^5 

- .00001 

-0.00001 

-0.00001 

-0.000008 

-0.00001 

TABLE  8 

METHODOLOGIES  OF  POTENTIAL-FLOW  PROGRAMS 


